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Abstract:  Seagrasses  are  rooted  flowering  marine  plants  that  provide  a  va¬ 
riety  of  ecosystem  services  to  the  coastal  areas  they  colonize.  Attenuation 
of  currents  and  waves  and  sediment  stabilization  are  often  listed  among 
these  services.  Although  we  have  a  reasonably  good  understanding  of  how 
currents  affect  seagrasses  and  vice-versa,  less  is  known  about  interactions 
between  waves  and  seagrasses,  and  standard  methods  for  research  on 
waves  in  seagrass  systems  have  not  yet  been  established.  This  report  pre¬ 
sents  background  information  needed  to  inform  and  encourage  further 
studies  on  waves  in  seagrass  systems  from  both  field  and  modeling  per¬ 
spectives.  It  reviews  current  knowledge  of  waves  in  seagrass  systems,  en¬ 
compassing  field  and  laboratory  data  as  well  as  modeling  efforts.  It  then 
describes  various  methods  for  measuring  waves  in  seagrass  colonized  ar¬ 
eas  and  modeling  the  dynamics  of  wave-seagrass  interactions.  Standardi¬ 
zation  of  experimental  designs,  instrumentation,  analyses,  and  modeling 
approaches  to  allow  for  ready  comparison  between  studies  is  encouraged. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Citation  of 
trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products.  All  product  names 
and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to  be  construed  as  an  offi¬ 
cial  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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and  forth  with  waves.  Therefore,  the  depicted  canopy  height  is  actually  the  maximum 
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canopy  (B).  Note  that  the  meadow  only  occupies  a  small  fraction  of  the  water  column 
while  the  canopy  occupies  the  entire  water  column.  In  this  case,  the  canopy  is  formed  by 
the  reproductive  structures  of  a  seagrass  but  often  seagrass  leaves  can  occupy  the  entire 


water  column . 5 
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in  the  Pacific  Ocean,  off  Baja  California,  Mexico . 6 
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power,  unless  dissipation  is  significant.  The  wave  steepens  as  a  result  of  both  changes. 


and  eventually  breaks . 8 
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apparent  in  the  shadow  of  the  breakwater  (right),  where  crossing  of  waves  from  different 


directions  also  leads  to  interference  patterns . 11 

Figure  7.  Wave  attenuation  in  a  dense  (1,270  ±  92  shoots/m^)  Ruppia  maritima  bed  off 

Bishop’s  Flead  Point  in  Chesapeake  Bay,  MD.  Note  that  wave  attenuation  is  inversely 

related  to  water/tidal  level . 13 

Figure  8.  An  example  of  how  sampling  rate  affects  the  results  of  wave  sampling 
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(0).  In  contrast,  if  the  same  short  waves  (E)  are  sampled  at  a  higher  rate,  their 


characteristics  are  much  better  represented  by  the  data  collected  (F) . 26 
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Figure  10.  Visual  representation  of  the  wave  breakdown  process  undertaken  in  the  Fast 
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(Figure  11) . 
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Figure  11.  Nondirectional  wave  spectra  for  two  wave  bursts  recorded  on  the  eastern  side 
of  the  Virginia  portion  of  Chesapeake  Bay  between  Nassawadox  Creek  and  Hungars 
Creek.  The  data  were  coiiected  with  a  5-MHz  Acoustic  Doppier  Veiocimeter  (ADV), 
equipped  with  a  pressure  gauge,  and  mounted  to  a  bottom  ianding  tripod  depioyed  in  3-m 
water  depth.  Burst  number  475  (biue  iine)  shows  a  sharp  peak  at  0.2  Hz  with  a  sharp 
deciine  in  energy  at  higher  frequencies,  representing  a  narrower,  weii-defined  wave  train. 

Burst  number  599  (red  iine)  shows  a  flatter  response  with  a  broad  peak  and  iess  faii-off 
of  energy  at  higher  frequencies,  representing  a  more  chaotic,  choppy  sea  state . 31 
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iegend  of  Figure  11 . 33 
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the  same  depth.  The  right  panei  gives  an  exampie  of  how  wave  sensors  can  be  depioyed 

at  the  same  depth.  Red  ovais  represent  essentiai  wave  gauges  whiie  yeiiow  stars 

represent  suggested  iocations  for  additionai  sensors . 34 
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a  month.  The  shoot  density  of  Ruppia  maritima  bed  is  about  1536, 1270,  and  1968/m2 
for  May  (a),  June  (b),  and  October  (c),  respectiveiy,  and  the  corresponding  canopy  height  is 
0.4, 1.0,  and  0.4  m.  The  modei  is  caiibrated  by  use  of  October  data  (see  text).  The  1:1 
iine  represents  no  wave  attenuation.  From  the  iinear  regression  of  fieid  data,  wave 
attenuation  is  higher  in  June  than  in  May  and  October  due  to  the  presence  of  iong 
reproductive  shoots . 
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1  Introduction 

Seagrasses  and  their  fluid  environment 

Seagrasses  are  marine  flowering  plants.  They  differ  from  macroalgae  (also 
known  as  seaweed)  as  these  do  not  flower  but  reproduce  via  spores.  Sea¬ 
grasses  also  have  a  root  and  true  rhizome  system  while  macroalgae  only 
have  rhizoids.  As  a  result,  macroalgae  tend  to  attach  to  hard  substrates 
such  as  shell  and  rocks  while  most  seagrasses  colonize  soft  substrates  such 
as  sand.  However,  some  seagrasses,  such  as  those  of  the  genus  Phyl- 
lospadix  and  sometimes  Posidonia  attach  to  rocks  and  are  exposed  to  rela¬ 
tively  high  wave  energy.  In  contrast,  seagrasses  colonizing  soft  substrates 
are  usually  found  in  more  quiescent  (i.e.,  wave-sheltered)  areas. 

Seagrasses  are  found  in  coastal  areas  around  the  world  except  Antarctica 
(Green  and  Short  2003)  and  are  limited  in  the  depths  and  areas  they  colo¬ 
nize  by  light  availability  (Dennison  et  al.  1993)  and  physical  and  geological 
parameters  (Koch  2001),  respectively.  Seagrasses  require  relatively  high 
light  levels  (11  to  20  percent  of  surface  light)  (Dennison  et  al.  1993;  Duarte 
1991;  Gallegos  and  Kenworthy  1996;  Kenworthy  and  Fonseca  1996)  when 
compared  to  other  marine  photosynthetic  organisms  such  as  phytoplank¬ 
ton,  which  only  requires  1  percent  surface  light.  As  a  result,  seagrasses  are 
generally  found  in  shallow  waters  (<  90  m)  (Duarte  1991)  also  influenced 
by  waves.  This  can  limit  the  geographical  distribution  of  seagrasses  as 
most  species  require  relatively  sheltered  conditions  to  thrive  (Fonseca  and 
Bell  1998;  Robbins  and  Bell  2000;  Koch  2001).  Excessively  sheltered  con¬ 
ditions  such  as  those  found  in  stagnant  ponds  or  shoreward  of  breakwaters 
can  also  be  detrimental  to  seagrasses  (Koch  2001)  and  other  organisms 
(Martin  et  al.  2005).  An  optimum  balance  between  light  availability  and 
moderate  wave  climate  is  often  found  in  estuaries,  embayments,  and 
coastal  lagoons.  Seagrasses  are  quite  abundant  in  these  areas.  Adjacent 
animal  communities  and/or  bathymetric  features  can  also  promote  sea- 
grass  establishment  via  wave  attenuation.  For  example,  extensive  seagrass 
beds  are  often  found  shoreward  of  coral  reefs  and  sandbars  where  wave 
conditions  are  relatively  quiescent  (see  review  by  Fonseca  1996). 

The  importance  of  seagrasses  in  coastal  systems  has  been  recognized  since 
early  settlers  used  their  seeds  as  a  food  source  and  their  leaves  as  insula¬ 
tion  material  (Wyllie-Echeverria  et  al.  2002).  Although  seagrasses  are  no 


ERDCTR-06-15 


2 


longer  used  for  these  purposes,  their  capacity  to  remove  nutrients  from  the 
water  column  has  recently  been  valued  $19,000  /hectare/  year  (Costanza 
et  al.  1997).  Other  ecosystem  services  provided  by  seagrasses  include  serv¬ 
ing  as  habitat  for  a  variety  of  ecologically  and  economically  valuable  spe¬ 
cies,  providing  food  to  associated  organisms  directly  (to  turtles)  and  indi¬ 
rectly  (via  detritus  and  epiphytes),  stabilizing  the  sediments  they  colonize, 
and  sequestering  carbon  (Green  and  Short  2003). 

Interactions  between  seagrasses  and  water  flow 

In  their  natural  environment,  seagrasses  are  exposed  to  wind-driven  cur¬ 
rents,  tides,  waves  and  wave-driven  currents  (Koch  2001).  While  these 
hydrodynamic  processes  affect  seagrasses,  seagrasses  also  affect  these  hy¬ 
drodynamic  processes.  A  feedback  mechanism  appears  to  develop  between 
seagrasses  and  their  fluid  environment  (Koch  et  al.  2006).  For  example,  as 
a  newly  established  seagrass  bed  attenuates  currents  and  waves,  fine  and 
organic  particles  tend  to  be  trapped  increasing  the  nutrient  availability  to 
the  plants  (Kenworthy  et  al.  1982).  As  a  result,  seagrasses  can  grow  faster 
and  more  robust  (Short  1987).  However,  excessive  attenuation  of  currents 
and  waves  by  extremely  dense  seagrass  beds  or  geomorphological  features 
may  lead  to  sediment  organic  contents  so  high  that  they  become  detrimen¬ 
tal  to  seagrasses.  Organic  contents  higher  than  5  percent  have  been  sug¬ 
gested  as  detrimental  to  seagrasses  (Koch  2001;  Kemp  et  al.  2004).  Addi¬ 
tionally,  excessively  weak  currents  and  waves  can  also  lead  to  limiting  leaf 
diffusive  boundary  layer  conditions  (Fonseca  and  Kenworthy  1987).  Cur¬ 
rents  below  5  cm/sec  have  been  suggested  as  limiting  to  seagrasses  due  to 
a  reduction  in  the  flux  of  carbon  and  nutrients  to  the  leaf  surface  (Koch 
1994)-  In  contrast,  in  areas  with  high  wave  exposure  and  strong  currents, 
seagrass  may  be  damaged  due  to  excessive  sediment  transport,  which  does 
not  allow  seeds  to  become  established,  or  eroding/burying  existing  sea¬ 
grass  beds.  As  a  result,  wave-  or  current-exposed  areas  tend  to  have  patchy 
seagrasses  or  are  unvegetated  (Fonseca  and  Bell  1998;  Hovel  et  al.  2002; 
Krause-Jensen  et  al.  2003),  and  an  intermediate  flow  regime  may  be  op¬ 
timal  for  seagrass  growth  and  development  (Koch  et  al.  2006). 

One  of  the  most  complex  water  flow  systems  is  that  in  vegetated  areas 
(Raupach  et  al.  1991,  Finnigan  2000).  Not  only  can  water  flow  affect  sea¬ 
grasses  and  seagrasses  affect  water  flow  (as  previously  described)  but  sea¬ 
grasses  and  water  flow  may  interact  in  highly  coupled,  nonlinear  ways. 
Kelvin-Helmholtz  instabilities  that  generate  large  coherent  vortices  at  the 
interface  between  the  canopy  and  the  overlaying  water  column  lead  to 
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wave-like  oscillations  referred  to  as  monami  (Ackerman  and  Okubo  1993). 
The  vortices  and  the  elastic  and  buoyant  seagrass  leaves  interact  causing 
the  canopy  to  wave  in  a  coherent  manner.  These  vortices  can  penetrate 
into  the  canopy  enhancing  vertical  transport  between  the  water  column 
and  the  canopy  (Ghisalberti  and  Nepf  2002).  As  a  result,  recruitment  of 
organisms  (Grizzle  et  al.  1996),  sediment  dynamics  and  nutrient  acquisi¬ 
tion  by  the  plants  can  be  enhanced. 

Due  to  the  complexities  of  water  flow  in  vegetated  systems,  hydrodynam¬ 
ics  of  seagrasses  is  a  relatively  unstudied  field  (Koch  2001).  Although  it  is 
beginning  to  be  understood  how  currents  move  through  seagrass  beds  and 
how  in  turn  seagrasses  affect  and  are  affected  by  currents  (Koch  et  al. 
2006),  understanding  about  waves  in  seagrass  systems  is  in  its  infancy 
(Fonseca  and  Cahalan  1992;  Fonseca  1998). 

The  purpose  of  this  publication  is  to  provide  engineers  and  scientists  with 
the  background  information  needed  to  encourage  further  studies  of  waves 
in  seagrass  systems  from  a  field  as  well  as  a  modeling  perspective.  Signifi¬ 
cant  progress  still  needs  to  be  made  regarding  how  waves  propagate 
through  seagrass  beds,  how  seagrasses  stabilize  the  sediments  they  colo¬ 
nize  and  possibly  also  the  adjacent  shoreline,  how  much  wave  energy  is 
tolerated  by  different  seagrass  species,  and  how  much  wave  energy  is 
beneficial  for  seagrass  systems,  just  to  cite  a  few.  This  report  begins  with  a 
review  the  current  knowledge  of  waves  in  seagrass  systems  including  field 
and  laboratory  data  as  well  as  modeling  efforts.  Then  methods  are  de¬ 
scribed  that  can  be  used  to  further  study  waves  in  seagrass  colonized  ar¬ 
eas.  It  concludes  with  the  latest  models  being  developed  to  understand  the 
dynamics  between  waves  and  seagrasses. 
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2  Background 

Seagrass  terminology 

Although  seagrass  morphology  closely  resembles  that  of  terrestrial 
grasses,  botanically,  seagrasses  are  more  closely  related  to  lilies  than 
grasses.  Seagrasses  have  a  horizontal  rhizome  linking  clusters  of  leaves 
referred  to  as  shoots,  and  roots  are  usually  found  at  each  shoot  (Figure  i). 
Seagrasses  are  true  flowering  plants.  As  a  result,  they  also  have  flowers 
and  seeds  (Figure  i).  The  flowers  of  some  species  (e.g.,  Halodule  wrightii, 
Thalassia  testudinum)  are  found  near  the  sediment  surface  while  other 
seagrass  species  (e.g.,  Zostera  marina,  Ruppia  maritima),  when  repro¬ 
ductive,  form  long  vertical  stems  that  can  occupy  most  of  or  the  entire  wa¬ 
ter  column.  When  most  of  the  plant  biomass,  independent  if  in  the  form  of 
reproductive  or  vegetative  shoots,  occupies  a  large  portion  of  the  water 
column,  the  vegetation  is  often  referred  to  as  a  canopy  (Figure  2).  In  con¬ 
trast,  when  most  of  the  seagrass  biomass  is  found  near  the  bottom,  they 
are  often  referred  to  as  meadows  (Figure  2). 


Figure  1.  Nomenclature  commonly  used  to  describe  parts  of  seagrasses  and  attributes  of  the 
canopy  they  form.  Note  that  the  leaves  bend  in  the  direction  of  currents  and  flap  back  and 
forth  with  waves.  Therefore,  the  depicted  canopy  height  is  actually  the  maximum  canopy 

height. 


ERDCTR-06-15 


5 


Figure  2.  Diagram  depicting  the  difference  between  a  seagrass  meadow  (A)  and  seagrass 
canopy  (B).  Note  that  the  meadow  oniy  occupies  a  smaii  fraction  of  the  water  coiumn  whiie 
the  canopy  occupies  the  entire  water  coiumn.  in  this  case,  the  canopy  is  formed  by  the 
reproductive  structures  of  a  seagrass  but  often  seagrass  ieaves  can  occupy  the  entire  water 

coiumn. 

When  quantifying  the  amount  of  seagrass  present  in  an  area,  quite  often 
the  number  of  shoots  in  an  area  (e.g.,  25  x  25  cm)  is  counted  and  ex¬ 
pressed  per  m2  (i.e.,  shoots/m^).  This  leads  to  a  parameter  called  “shoot 
density”.  This  value  varies  between  species.  For  example,  a  dense  Ruppia 
maritima  bed  may  have  more  than  3,000  shoots/m^  while  a  dense  Zostera 
marina  bed  may  only  have  2,000  shoots/m^.  This  appears  to  be  due  to  the 
biomass  occupied  by  the  leaves  or  leaf  area  (note  that  Zostera  has  wider 
leaves  than  Ruppia).  Canopy  height  is  also  an  important  seagrass  parame¬ 
ter  when  evaluating  the  impact  of  seagrasses  on  water  flow.  This  parame¬ 
ter  is  usually  obtained  by  averaging  the  tallest  two-thirds  of  the  leaves 
(Figure  1).  The  smallest  seagrasses  (genus  Halophila)  are  only  2  or  3  cm 
tall  while  the  largest  seagrasses  in  the  United  States  can  reach  2  m  in 
length  (nonreproductive  Zostera  marina  and  Phyllospadix  sp.).  When  re¬ 
productive,  other  seagrasses  such  as  Ruppia  maritima  can  also  reach  1  m 
in  length.  Although  canopy  height  measured  in  the  form  of  the  tallest 
leaves  or  reproductive  structure  is  easy  to  quantify,  it  may  not  be  a  good 
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representation  of  the  true  canopy  height  observed  in  situ.  When  exposed 
to  currents  and  waves,  seagrass  leaves  tend  to  bend  in  the  direction  of  the 
flow  (Figure  3)  leading  to  a  canopy  height  smaller  than  the  leaf  or  repro¬ 
ductive  shoot  length.  Fonseca  et  al.  (1982)  described  the  angle  of  bending 
of  Zostera  marina  shoots  as  a  function  of  current  velocity  in  a  flume  by  the 
following  equation: 


Y  =  62.8e-““'’“ 


(1) 


where  Y  is  the  bending  angle  measured  from  the  horizontal  and  u  is  the 
current  velocity  in  cm  s  ^.  Fonseca  and  Fisher  (1986)  also  describe  seagrass 
canopy  compression  (canopy  height/water  depth)  as  a  function  of  Darcy’s 
friction  factor.  Due  to  the  complexity  of  modeling  flexible  structures  in  the 
aquatic  environment,  stiff  cylinders  are  commonly  considered  in  models 
(see  Nepf  1999)  which  assume  that  canopy  height  equals  leaf  length.  As 
seen  from  the  Equation  1,  this  is  an  overestimation  of  canopy  height  and 
its  effects  on  water  flow. 


Figure  3.  Phyllospadix  swaying  back  and  forth  with  the  waves  at  a  shaiiow  site  (2  to  3  m)  in 
the  Pacific  Ocean,  off  Baja  Caiifornia,  Mexico. 

Nearshore  waves:  A  brief  summary  of  important  characteristics 

In  the  most  general  sense,  a  wave  may  be  defined  as  a  disturbance  that 
travels  in  at  least  one  spatial  dimension  through  time,  and  to  a  good  ap¬ 
proximation  transports  energy  without  transporting  mass.  Waves  that 
travel  on  the  surface  of  the  water  are  known  as  surface  gravity  waves;  these 
are  the  waves  that  most  directly  affect  nearshore  seagrass  beds.  The  waves 
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that  are  the  main  focus  of  this  review  are  those  that  are  generated  (some¬ 
where)  by  the  wind  and  travel  in  the  direction  of  the  wind.  Because  the 
time  and  space  scales  of  generation  of  these  waves  are  long,  they  may  be 
approximated  mathematically  by  a  continuous  sine  (or  cosine)  function,  so 
the  surface  elevation  p  of  a  wave  traveling  over  a  distance  x  through  time  t 
can  be  expressed  as: 


.  .  H  .  .2nx  2nt. 

ri{x,t)  =  —sm{— - — )  (2) 

^  1 

In  Equation  2,  H  is  defined  as  the  local  wave  height,  the  vertical  distance 
between  a  crest  (the  highest  point  of  the  wave)  and  the  adjacent  trough 
(the  lowest  point  of  the  wave);  Figure  4.  L  is  defined  as  the  wavelength  of 
the  wave,  the  distance  between  successive  crests  in  the  direction  of  wave 
travel.  T  is  defined  as  the  period  of  the  wave,  the  time  between  the  pas¬ 
sage  of  two  successive  crests  measured  by  a  fixed  observer.  The  phase 
speed  (the  speed  at  which  a  wave  crest  travels)  is  defined  by  c=L/T.  Equa¬ 
tion  2  is  known  as  the  linear  or  Airy  wave  approximation,  and  it  is  quite 
approximate  for  a  wind  sea  as  it  approaches  the  shoreline.  Real  waves  in¬ 
clude  variations  in  wave  height,  period,  and  direction,  and  deviate  from 
the  sinusoidal  shape  in  shallow  water.  However,  it  is  remarkable  how  well 
this  expression  and  the  relationships  that  result  from  it  describe  basic 
wave  dynamics  in  the  nearshore.  Nearshore  waves  may  also  be  generated 
by  transient  disturbances  such  as  boating  traffic  (boat  waves),  or  by  inter¬ 
actions  between  breaking  waves  and  the  shoreline  geometry  (infragravity 
waves),  but  these  wave  types  are  not  well  described  by  Equation  2  so  are 
not  considered  further  here. 

The  most  important  aspect  of  wave  dynamics  is  the  transport  of  energy. 
The  energy  E  of  a  wave  averaged  over  a  wave  period  is  defined  by: 

E  =  \pgH^  (3) 

O 

where  p  is  the  density  of  water  and  g  is  the  gravitational  constant,  in  units 
of  Joules  (or  Watt-sec)  per  unit  of  surface  area.  This  energy  is  transported 
in  the  direction  of  wave  propagation  at  the  group  velocity  Cg  of  the  wave, 
such  that  the  rate  of  wave  energy  transport,  also  called  the  wave  energy 
flux  or  wave  power,  is  given  by: 
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Ef=c^E 


in  watts  per  unit  of  wave  crest  length. 


(4) 


Figure  4.  Schematic  diagram  of  a  shoaiing  wave  as  it  approaches  shore.  Variabie  names  are 
defined  in  the  text.  As  the  depth  h  becomes  shaiiower  nearshore,  the  speed  of  the  wave 
decreases.  The  wave  period  cannot  change,  so  in  order  for  the  speed  to  decrease  the 
waveiength  must  decrease.  The  wave  height  must  aiso  increase  to  conserve  wave  power, 
uniess  dissipation  is  significant.  The  wave  steepens  as  a  result  of  both  changes,  and 

eventuaiiy  breaks. 

Note  that  Cg  is  not  necessarily  the  same  as  the  wave  phase  speed  c  shown 
in  Figure  4.  In  general,  the  group  velocity  is  given  by: 


1 


1  +  - 
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(5) 


where  L  is  given  by  the  solution  of  the  transcendental  function: 


Itt 

kTj 
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tanh 


2nh 

L 


(6) 


For  shallow-water  waves  with  h/L  <  0.04,  where  h  is  the  water  depth,  the 
phase  speed  and  group  velocity  are  the  same  and  depend  only  on  water 
depth: 


Cg=c  =  4^ 


(7) 
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For  deep  water  waves  with  h/L  > 


That  is,  the  energy  only  travels  half  as  fast  as  the  individual  waves  and 
both  depend  only  on  the  wave  period.  Between  these  limits,  for  transi¬ 
tional  or  intermediate  waves,  the  expressions  for  c  and  Cg  are  complex 
functions  of  both  depth  and  period  (see  e.g..  Dean  and  Dalrymple  1991). 

Figure  5  shows  the  group  velocity  Cg  for  wave  periods  and  depths  relevant 
for  most  nearshore  seagrass  beds.  Clearly,  the  behavior  of  waves  interact¬ 
ing  with  seagrass  beds  can  vary  significantly  depending  on  the  wave  peri¬ 
ods  and  depths  of  interest.  Note  that  for  long-period  (e.g.,  T  =  12  sec), 
shallow-water  waves  Cg  changes  rapidly  as  shore  is  approached.  These 
waves  are  strongly  affected  by  the  local  depth,  produce  strong  near-bottom 
wave  action,  and  are  subject  to  significant  steepening  (through  wave  shoal¬ 
ing)  and  breaking.  For  short-period  (e.g.,  T  =  2  sec)  waves  Cg  does  not 
change  much  until  the  waves  are  in  shallow  water.  These  waves  are  only 
weakly  affected  by  the  local  depth,  produce  weaker  near-bottom  wave  ac¬ 
tion  (for  equivalent  wave  height),  and  experience  much  less  steepening 
and  breaking  than  the  long-period  waves  in  the  same  water  depths.  Waves 
interact  with  the  bottom  based  on  the  relative  water  depth  (h/L),  so  long- 
period  waves  have  longer  wavelengths  and  interact  with  the  bottom  at 
much  deeper  depths. 

The  waves  of  interest  here  are  generated  by  the  action  of  the  wind  stress 
on  the  surface  of  the  water.  Waves  become  longer  and  higher  as  the  fetch 
and/or  the  duration  increase.  Fetch  is  the  over- water  distance  in  the  di¬ 
rection  of  the  wind,  and  duration  is  the  amount  of  time  the  wind  has  been 
blowing  from  that  direction.  Waves  generated  in  shallow  water  tend  to  be 
shorter  and  smaller  than  waves  generated  in  deep  water.  Swells  are  long 
waves  generated  in  deep  water  that  outrun  their  source  region,  sometimes 
traveling  thousands  of  kilometers  before  encountering  land.  They  tend  to 
be  unidirectional  (from  a  single  direction),  narrow  banded  (energy  is  con¬ 
centrated  in  a  narrow  range  of  periods),  and  to  travel  in  groups  of  larger 
waves.  Waves  generated  by  local  winds,  often  referred  to  as  chop  or  wind- 
sea,  are  more  chaotic  short  waves  with  a  broad  spectrum  of  wave  heights, 
periods,  and  directions  centered  on  mean  values  that  can  change  rapidly 
as  the  wind  changes.  Long-period  waves  (swells  and  seas)  tend  to  domi- 
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nate  nearshore  environments  with  open  ocean  exposure,  but  short-period 
wind  waves  usually  dominant  in  closed  or  semi-enclosed  water  bodies. 
Wave  power,  once  generated,  is  either  conserved  or  dissipated  as  a  wave 
approaches  shore  and  is  transformed  by  interactions  with  the  bottom 
slope,  seagrass  beds,  or  underwater  obstacles  such  as  breakwaters  or  off¬ 
shore  bars.  At  the  shoreline,  any  remaining  wave  energy  is  either  reflected 
back  or  dissipated  by  breaking.  Violent  dissipation  of  wave  energy  in  the 
nearshore  through  breaking  is  responsible  for  the  high  rates  of  sediment 
transport  and/or  shoreline  erosion  that  commonly  occur  along  wave- 
exposed  coasts. 


Water  depth  (m) 


Figure  5.  Group  velocities  predicted  by  linear  wave  theory  for  7=2-12  sec. 


Wave  transformation  patterns  are  often  complex  near  the  shore,  especially 
in  the  presence  of  complex  bathymetry  and/or  shoreline  geometry,  or  sig¬ 
nificant  local  changes  in  wave  drag  (Figure  6).  Waves  approaching  the 
shoreline  at  an  angle  will  bend  toward  the  shoreline  as  the  group  velocity 
nearest  the  shore  slows  down;  this  process  is  known  as  refraction.  Waves 
encountering  an  obstacle  of  limited  length  (a  breakwater  or  seawall,  for 
example)  will  be  entirely  or  partially  blocked  by  the  obstacle,  but  wave  en¬ 
ergy  will  spread  into  the  shadow  of  the  obstacle  from  beyond  its  edges  via 
diffraction.  A  marked  local  change  in  bottom  roughness/drag,  such  as 
produced  by  an  oyster  reef  or  a  seagrass  bed,  will  have  some  of  the  same 
effects  as  a  breakwater  except  that  complete  blocking  of  wave  energy  is  re¬ 
placed  by  more  gradual  energy  dissipation.  This  allows  more  of  the  inci¬ 
dent  wave  energy  to  continue  propagating  past  the  reef/bed  than  for  a 
breakwater.  In  all  cases,  changes  in  wave  height  are  determined  by  the  re¬ 
quirement  that  the  wave  energy  flux  budget  be  balanced,  such  that  the  to- 
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tal  energy  flux  between  two  wave  rays  (lines  perpendicular  to  the  crests  in 
the  direction  of  propagation)  is  either  conserved  or  dissipated.  The  com¬ 
bination  of  all  of  these  wave  transformation  effects  as  waves  transition 
from  deep  to  shallow  water  virtually  mandates  the  use  of  a  numerical 
model  for  anything  more  than  the  simplest  of  problems. 


Figure  6.  Illustrations  of  refraction  (left  panel)  and  diffraction  (right  panel).  The  solid  lines 
represent  the  wave  crests,  the  short  dashed  lines  represent  depth  contours,  and  the  long 
dashed  lines  represent  wave  rays.  Turning  and  spreading  of  wave  rays  are  apparent  in  both 
cases,  resulting  in  a  redirection  and  decrease  of  wave  power.  This  decrease  is  most  apparent 
in  the  shadow  of  the  breakwater  (right),  where  crossing  of  waves  from  different  directions  also 

leads  to  interference  patterns. 
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3  How  Seagrasses  Affect  Waves  and  Wave- 
Dependent  Processes 

Wave  attenuation 

Seagrasses  tend  to  attenuate  waves,  especially  in  shallow  waters  where 
waves  interact  with  the  seagrass  canopy  the  same  way  as  shallow  waves 
interact  with  the  bottom.  Even  so,  most  waves  in  seagrass  beds  seem  to  be 
intermediary  in  nature  (i.e.,  0.04  <  h/L  <  0.5).  The  degree  of  wave  at¬ 
tenuation  observed  in  a  seagrass  bed  is  directly  related  to  the  fraction  of 
the  water  column  occupied  by  the  vegetation.  Wave  attenuation  is  highest 
when  seagrasses  occupy  a  large  portion  (>50  percent)  of  the  water  column 
(Ward  et  al.  1984;  Fonseca  and  Cahalan  1992;  Figure  7).  A  flume  study 
measured  wave  attenuations  between  20  and  76  percent  over  i-m  length 
when  the  plants  were  occupying  the  entire  water  depth  (Fonseca  and  Ca¬ 
halan  1992),  whereas  field  studies  measured  values  between  1.6  and  80 
percent  (Koch  1996;  Prager  and  Halley  1999).  This  wide  variation  in  field 
data  appears  to  be  related  to  tidal  fluctuations,  i.e.  fraction  of  the  water 
column  occupied  by  the  vegetation:  at  high  tide,  wave  attenuation  is 
smaller  than  at  low  tide  (Figure  7).  Reduction  in  wave  energy  is  also  ob¬ 
served  in  relatively  deep  beds  exposed  to  long  waves.  For  example,  Ver- 
duin  and  Backhaus  (2000)  observed  that  15-sec  waves  were  attenuated  by 
anAmphibolis  antarctica  bed  at  5-m  depth  (note  that  h/L  =  0.33,  i.e.,  in¬ 
termediate  waves). 

Flooding  is  usually  associated  with  severe  storm  events.  Therefore,  it  could 
be  speculated  that  seagrass  beds  will  be  less  effective  in  attenuating  waves 
(canopy  occupies  a  smaller  fraction  of  the  water  column)  when  this  ecosys¬ 
tem  service  is  needed  the  most.  However,  storms  also  generate  longer 
waves.  Therefore,  theoretically,  seagrasses  may  still  attenuate  waves  dur¬ 
ing  storm  events  when  flooding  occurs.  This  complex  relationship  between 
storm  waves,  flooding,  and  seagrass-induced  wave  attenuation  still  needs 
to  be  verified. 
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Figure  7.  Wave  attenuation  in  a  dense  (1,270  ±  92  shoots/m^)  Ruppia  maritima  bed  off 
Bishop’s  Head  Point  in  Chesapeake  Bay,  MD.  Note  that  wave  attenuation  is  inversely  related 

to  water/tidal  level. 


Sediment  stabilization 

An  ecosystem  service  commonly  listed  for  seagrasses  is  sediment  and 
shoreline  stabilization.  Seagrasses  effectively  reduce  current  velocities  and 
attenuate  waves  (Fonseca  et  al.  1982;  Fonseca  and  Fisher  1986;  Gambi  et 
al.  1990;  Fonseca  1996;  Koch  1996;  Wallance  and  Cox  1997;  Koch  and 
Gust  1999;  Verduin  and  Backhaus  2000;  Peterson  et  al.  2004)  and,  as  a 
result,  sediment  erosion  is  decreased  and  particle  deposition  is  enhanced 
in  seagrass  beds.  Consequently,  sediments  in  vegetated  areas  are  usually 
finer  and  more  organic  than  in  adjacent  unvegetated  areas  (Kenworthy  et 
al.  1982;  Fonseca  and  Koehl  2006).  Quiescent  hydrodynamic  conditions 
are  likely  to  also  contribute  to  the  high  recruitment  of  larvae  and  juvenile 
stages  of  many  organisms  such  as  crustaceans  and  fishes. 

Due  to  their  capacity  to  alter  their  environment,  seagrasses  have  been  re¬ 
ferred  to  as  ecosystem  engineers  (Jones  et  al.  1994, 1997;  Thomas  et  al. 
2000).  Caution  is  needed  though  as  seagrasses  do  not  always  fit  this  label. 
In  highly  wave-exposed  sites  where  seagrasses  do  not  attenuate  water  flow 
as  effectively  as  in  unidirectional  (tidal)  flows  (Koch  and  Gust  1999), 
sediment  characteristics  within  and  outside  seagrass  beds  differed  little  to 
none  (Hoskin  1983;  Edgar  and  Shaw  1991;  Koch  1999;  van  Keulen  and 
Borowitzka  2001).  Actually,  in  some  cases,  sediment  in  a  vegetated  area 
can  be  coarser  than  in  the  adjacent  unvegetated  area  (Koch  1993;  Fonseca 
and  Koehl  2006;  Koch  et  al.  2006). 
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4  How  Waves  Affect  Seagrasses 

Negative  effects 

Erosion/dislodgement 

Landscape  patterns  observed  in  seagrass  habitats  are  often  associated  with 
hydrodynamic  disturbances  (Fonseca  et  al.  1983;  Fonseca  and  Bell  1998; 
Hemminga  and  Duarte  2000).  As  described  in  Koch  et  al.  (2006),  in  areas 
disturbed  by  high  wave  energy,  seagrass  ecosystems  can  be:  (1)  non¬ 
existent  (Dan  et  al.  1998);  (2)  depth  restricted  (when  sufficient  light  is 
available,  seagrasses  colonize  areas  below  the  maximum  wave  penetration 
depth  (Krause- Jensen  et  al.  2003;  Middelboe  et  al.  2003);  (3)  dominated 
by  more  robust  species  (e.g.,  Amphibolis  grifftthii  and  Posidonia  cori- 
acea);  or  (4)  patchier  as  the  disturbance  of  high  waves  may  hinder  the  lat¬ 
eral  expansion  of  some  seagrass  beds  (Kendrick  et  al.  2000;  Frederiksen  et 
al.  2004).  In  contrast,  in  sheltered  waters,  seagrass  meadows  tend  to  be 
more  continuous  and  are  colonized  by  relatively  more  fragile  species  (e.g., 
Posidonia  spp)  (Kirkman  and  Kuo  1990). 

Recent  developments  (Koch,  in  preparation)  suggest  that  seagrasses  are 
resilient  when  exposed  to  waves  (also  see  the  following  section).  Their 
massive  loss  occurs  when  sediment  becomes  unstable  due  to  erosion.  Sea¬ 
grasses  can  also  be  lost  due  to  excessive  sedimentation  (Marba  et  al.  1994; 
Walker  et  al.  1996;  Bridgwood  2002;  Paling  et  al.  2003,  van  Keulen  and 
Borowitzka  2003;  Frederiksen  et  al.  2004).  Sand  waves,  a  result  of  water 
flow,  can  lead  to  the  deposition  of  tens  of  centimeters  over  periods  of 
hours  (Paling  et  al.  2003).  The  degree  to  which  these  large  amounts  of 
sediment  negatively  affect  the  seagrasses  creating  unvegetated  patches  de¬ 
pends  on  their  tolerance  to  sedimentation,  the  amount  of  sediment  depos¬ 
ited,  and  the  period  the  plants  remain  buried.  Some  seagrasses  (e.g., 
Halodule  wrightii)  are  able  to  survive  as  long  as  the  sediment  is  removed 
in  a  matter  of  weeks  (Phillips  1980)  while  others  (e.g.,  Z.  marina)  seem  to 
have  little  or  no  tolerance  to  sedimentation  regardless  of  the  sediment  type 
(Mills  and  Fonseca  2003). 

Mechanical  tissue  damage 

In  wave  swept  environments,  a  long  flexible  shape  (see  Figure  3)  can  mini¬ 
mize  the  forces  exerted  on  the  roots  of  seagrasses.  The  long  leaves  tend  to 
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move  back  and  forth  without  fully  extending,  thereby  minimizing  the  drag 
imposed  on  the  roots  (Koehl  1984). 

The  field  of  seagrass  biomechanics  is  mostly  unstudied.  Kopp  (1999) 
showed  that  strength  of  vegetative  Z.  marina  shoots  varies  seasonally  with 
the  most  vulnerable  month  (highest  breakage)  in  September  when  the 
plants  begin  to  senesce  (i.e.,  when  plants  naturally  lose  their  leaves  due  to 
shorter  days  and/or  cooler  weather).  During  this  period,  more  than  50 
percent  of  the  leaves  are  unable  to  withstand  4.7  m/sec  (current  speed  in 
seagrass  beds  is  usually  around  0.1  m/sec  during  calm  conditions).  Patter¬ 
son  et  al.  (2001)  showed  that  reproductive  shoots  of  Z.  marina  are  as 
tough  as  macroalgae  and  that,  among  a  population,  a  few  reproductive 
shoots  are  extremely  tough.  These  may  persist  even  during  hurricane 
forces  ensuring  the  survival  of  at  least  some  individuals  (assuming  that 
sediment  erosion  will  not  compromise  the  entire  bed).  Fonseca  (1998) 
suggested  that  the  force  needed  to  damage  live  Z.  marina  may  be  achieved 
only  infrequently  under  extreme  storm  conditions. 

Turbidity 

Seagrasses  need  relatively  high  light  levels  in  order  to  thrive.  Therefore, 
anything  that  increases  water  turbidity  is  detrimental  to  seagrasses.  While 
sediments  tend  to  be  deposited  in  seagrass  beds  under  calm  conditions 
(see  “Sediment  stabilization”  section  in  Chapter  3),  some  of  these  particles 
may  be  resuspended  by  waves  during  storm  events.  Seagrasses  have  accli¬ 
mated  to  such  relatively  short  (hours)  pulsed  high  turbidity  events.  Tur¬ 
bidity  events  that  last  weeks  or  months  (e.g.,  shoreline  erosion,  excessive 
river  runoff,  phytoplankton  bloom)  are  the  ones  that  lead  to  the  loss  of 
seagrasses  (Moore  et  al.  1997).  Additionally,  wave-induced  turbidity  could 
potentially  be  detrimental  to  seagrasses  when  islands  are  lost  or  peninsu¬ 
las  breached  due  to  erosion.  This  leads  to  a  permanent  increase  in  wave 
energy  in  areas  that  were  relatively  protected  in  the  past.  As  a  result,  tur¬ 
bidity  is  expected  to  increase,  especially  in  the  beginning  when  fine  parti¬ 
cles  deposited  during  more  quiescent  conditions  are  resuspended. 

Positive  effects 

Flushing 

Waves  cause  seagrass  leaves  to  move  back  and  forth.  When  leaves  bend 
horizontally,  they  tend  to  isolate  the  sediment  and  water  within  the  canopy 
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from  the  water  column.  Just  a  fraction  of  a  second  later,  when  the  leaves 
are  moved  to  a  more  vertical  position  by  a  passing  wave,  the  canopy  opens 
and  the  exchange  between  the  water  column  and  the  seagrass  canopy  in¬ 
creases.  The  leaves  then  return  to  a  more  horizontal  position  starting  the 
cycle  over  again.  This  constant  leaf  motion  leads  to  extensive  flushing  of 
the  seagrass  bed  possibly  maximizing  the  flux  of  nutrients  and  carbon  to 
the  plant  surface  (Koch  and  Gust  1999). 

Reduced  self-shading  and  light-flecks 

Seagrasses  exposed  to  unidirectional  flow  (e.g.,  tidal  currents)  experience 
a  high  degree  of  self-shading  as  leaves  lay  on  top  of  each  other  for  ex¬ 
tended  periods  of  time.  In  contrast,  wave-induced  leaf  flapping  allows 
specks  of  light  to  penetrate  the  seagrass  canopy  at  the  frequency  of  flap¬ 
ping.  As  a  result,  higher  productivity  is  expected  in  wave-dominated  sys¬ 
tems.  Additionally,  wave  crests  focus  the  light  that  reaches  the  water  sur¬ 
face  leading  to  short  (fraction  of  a  second)  specks  of  high  light  levels 
referred  to  as  light-flecks  (Wing  and  Patterson  1993).  The  frequency  of 
light-flecks  resembles  that  of  the  passing  waves  (Koch  et  al.  2006).  When 
seagrasses  grow  in  shallow  waters,  they  may  benefit  from  these  light- 
flecks.  Although  their  effect  on  seagrass  productivity  was  not  yet  tested, 
phytoplankton  and  macroalgae  productivity  is  enhanced  by  light-flecks 
(Dromgoole  1988;  Greene  and  Gerard  1990;  Wing  and  Patterson  1993; 
Wing  et  al.  1993).  It  is  assumed  that  seagrasses  also  benefit  from  light- 
flecks. 

Pulses  in  the  diffusive  boundary  layer 

The  diffusive  boundary  layer  on  seagrass  leaves  serves  as  a  barrier  to  the 
free  flux  of  nutrients  and  carbon  to  the  plant  surface.  As  a  result,  relatively 
stagnant  conditions  reduce  the  productivity  of  seagrasses  (Koch  1994). 
Wave-induced  oscillatory  flows  tend  to  disrupt  the  diffusive  boundary 
layer  for  a  fraction  of  a  second.  This  process  supplies  pulses  of  molecules 
to  the  blade  surface  replenishing  nutrients  and  carbon  needed  for  plant 
growth  (Nikora  et  al.  2002).  Under  wave-dominated  conditions,  these 
pulses  occur  on  a  regular  basis  enhancing  productivity  (Stevens  and  Hurd 
1997)- 
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Reduced  epiphytic  cover 

Epiphytes  are  microalgae  that  grow  on  the  surface  of  seagrasses  blocking 
light  and,  when  in  excess,  leading  to  the  death  of  the  plants.  Epiphytic 
loading  usually  increases  with  nutrient  availability.  It  has  been  speculated 
that  waves  and  currents  can  lead  to  reduced  epiphytic  growth  on  seagrass 
leaves;  little  data  are  available  to  confirm  this  hypothesis.  Pinckney  and 
Micheli  (1998)  observed  no  difference  in  total  epiphyte  biomass  in  a  wave- 
exposed  and  a  sheltered  seagrass  habitat  but  noted  that  diatoms,  coralline, 
and  some  filamentous  algae  dominate  under  wave-exposed  conditions 
while  blue  green  and  other  filamentous  algae  dominate  under  calm  condi¬ 
tions. 

Seagrass  restoration  in  wave-dominated  systems 

Comparative  studies  of  seagrass  restoration  across  wave  energy  regimes 
are  few.  Individual  studies  show  limited  success  in  areas  with  high  wave 
energy  mainly  due  to  mechanical  removal  of  the  plants  and  extensive 
sediment  transport  (van  Katwijk  and  Hermus  2000;  Paling  et  al.  2003; 
Campbell  and  Paling  2003).  Sediment  erosion  and  transport  can  be  re¬ 
duced  by  using  shell  (van  Katwijk  and  Hermus  2000)  or  mats  (when  these 
stay  in  place;  Campbell  and  Paling  2003)  thereby  increasing  restoration 
success.  Restoration  success  in  areas  of  high  wave  exposure  can  also  be  in¬ 
creased  when  using  larger  size  seagrass  plugs  (10  to  15  cm  in  diameter;  van 
Keulen  et  al.  2003)  although  some  seagrass  habitats  have  such  high  wave 
energy  (up  to  50  cm  S'^  orbital  velocities  near  the  sediment  surface)  that 
not  even  the  largest  seagrass  sod  units  (0.5  x  0.5  x  0.35  m)  remain  in  place 
(Paling  et  al.  2003).  Increased  turbidity  in  wave  exposed  areas  could  pos¬ 
sibly  also  limit  seagrass  restoration,  but  this  effect  can  potentially  be  me¬ 
diated  by  oysters  which  filter  suspended  particles  thereby  increasing  light 
availability  (Newell  and  Koch  2004). 

Except  for  a  few  studies,  including  the  ones  previously  mentioned,  com¬ 
parative  analyses  of  wave  effects  on  restoration  attempts  remain  largely 
anectodal  or  based  on  modeling  exercises  (Kelly  et  al.  2001).  Eonseca  et  al. 
(1998)  provide  guidance  on  selection  of  sites  and  planting  arrangement 
based  on  relative  exposure  of  sites  to  waves,  but  this  too  is  the  result  not  of 
manipulative  experiments  but  experience  and  correlations  derived  from 
study  of  natural  seagrass  landscapes  (Bell  et  al.  1997).  In  general,  if  sea¬ 
grass  patches  are  expected  to  form,  rather  than  continuous  cover,  then  the 
wave  climate  exceeds  the  requirement  for  restoration.  In  this  case,  fre- 
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quent  replanting  efforts  are  likely  to  be  necessary  and  may  require  as 
many  as  50  percent  of  the  original  plant  number. 
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5  Measuring  Waves  in  Seagrass  Systems 

Introduction 

Although  seagrass  habitats  are  influenced  by  currents  and  waves,  studies 
focusing  on  seagrass  hydrodynamics  are  few  (Koch  2001).  Those  that  ad¬ 
dress  this  topic  focus  mainly  on  currents.  The  study  of  waves  in  seagrass 
habitats  is  still  in  its  infancy  (Koch  et  al.  2006).  As  a  result,  methods  used 
to  quantify  the  wave  climate  and  wave  exposure  in  seagrass  beds  are  few. 
Valuable  suggestions  can  be  drawn  from  the  kelp  literature  as  the  study  of 
waves  in  these  systems  is  more  advanced  than  in  seagrass  beds.  Caution  is 
needed  though  as  waves  in  kelp  beds  are  usually  oceanic  (long)  in  nature 
while  those  in  seagrass  beds  are  usually  shorter.  Therefore,  waves  need  to 
be  recorded  at  a  higher  frequency  in  seagrass  habitats  than  in  kelp  beds. 

Cartographic  methods 

In  order  for  (most)  waves  to  develop,  wind  needs  to  blow  over  the  water. 
The  longer  the  distance  over  which  the  wind  can  blow  without  encounter¬ 
ing  obstructions,  i.e.,  the  longer  the  fetch,  the  larger  the  waves  can  de¬ 
velop.  As  a  result,  the  quantification  of  fetch  is  a  simple  method  of  ap¬ 
proximating  the  wave  exposure  at  a  site.  Cartographic  methods  are  based 
on  this  principle. 

Baardseth  index  of  wave  exposure 

This  index  has  not  yet  been  used  for  seagrasses  but  is  often  applied  in 
studies  involving  macroalgae  (see  Sjotun  and  Fredriksen  1995;  Ruuskanen 
et  al.  1999;  and  references  therein).  The  Baardseth  index  can  be  deter¬ 
mined  by  marking  the  study  site  on  a  nautical  chart  and  drawing  a  circle 
around  the  study  site.  The  circle  is  then  subdivided  into  40  sectors 
(9  deg  each)  and  the  number  of  sectors  without  any  land  obstructions  such 
as  peninsulas  and  island  is  counted.  A  high  value  indicates  a  high  degree  of 
wave  exposure  (i.e.,  few  land  obstructions  to  protect  the  site)  while  a  value 
of  o  represents  complete  protection  (i.e.,  site  protected  by  land  from  all 
sides).  The  radius  of  the  circle  depends  on  the  study  site  and  the  bottom 
topography.  For  example,  Ruuskanen  at  al.  (1999)  used  a  radius  of  7.5  km 
at  a  site  surrounded  by  small  islands.  In  contrast,  Sjotun  and  Fredriksen 
(1995)  used  a  series  of  circles/radiuses  (0.5,  7  and  100  km)  and  calculated 
the  Baardseth  index  for  each  radius  (ei,  62,  63)  in  an  area  surrounded  by 
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islands  but  also  affected  by  the  local  topography  (a  fjord).  Total  wave  ex¬ 
posure  (E)  was  then  estimated  using  the  following  equation: 

E  =  [ei +(10^2) +  (100^3)]. 100  (9) 

Although  simple  to  estimate,  the  Baardseth  index  does  not  take  into  ac¬ 
count  the  dominant  wind  direction  and  intensity  or  the  bottom  bathym¬ 
etry.  This  may  lead  to  an  overestimation  of  wave  exposure  especially  in  ar¬ 
eas  in  which  sills  and  sandbars  occur.  Sjotun  and  Fredriksen  (1995) 
considered  wind  in  the  Baardseth  index  by  assigning  a  relative  wind  force 
value  (e)  to  each  sector  based  on  mean  wind  force  and  frequency  over  the 
last  5  years: 
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where  Sn  is  the  number  of  unobstructed  sectors  in  a  given  direction,  Un  is 
the  number  of  observations  of  wind  from  a  given  direction  and  Fn  is  the 
average  strength  of  wind  from  a  given  direction. 

Relative  wave  exposure  indices 

Keddy  (1982)  took  into  account  that  a  body  of  water  affects  its  wave  gen¬ 
eration  (e.g.,  although  the  fetch  straight  up  a  river  may  be  long,  the  adja¬ 
cent  shorelines  limit  wave  development)  and  considered  effective  fetch  in 
the  estimation  of  wave  exposure  in  lakes  colonized  by  freshwater  plants. 
The  effective  fetch  is  estimated  by  measuring  fetch  along  four  lines  radiat¬ 
ing  out  from  either  side  of  a  compass  heading  at  increments  of  11.25  deg. 
Keddy’s  relative  wave  exposure  index  (REI)  also  takes  into  account  wind 
velocity,  the  dominant  wind  direction  and  exceedance  winds: 


8 

REI  =  Y^(Vi-xPi-XFi) 

i=i 


(11) 


where  i  =  compass  headings  (1-8  (N,  NE,  E,  etc.),  in  45-deg  increments), 
V  =  average  monthly  maximum  wind  speed  in  m/sec  (or  any  other  wind 
data  of  choice),  P  =  frequency  as  a  percentage  of  time  wind  occurred  from 
the  2th  direction,  and  F  =  effective  fetch  (m).  A  detailed  description  of  this 
index  and  how  it  relates  to  the  distribution  of  submersed  aquatic  vegeta¬ 
tion  in  lakes  can  be  found  in  Keddy  (1982).  It  appears  that  Keddy’s  wave 
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exposure  index  is  capable  of  predicting  submersed  plant  distribution  in 
areas  with  steep  slopes  (e.g.,  some  lakes)  where  waves  are  only  affected  by 
the  bottom  close  to  shore.  In  areas  dominated  by  extensive  shallow  flats 
this  index  does  not  appear  to  hold  because  waves  interact  with  the  bottom 
(i.e.,  shallow  waves)  and,  therefore,  this  leads  to  a  disconnection  between 
fetch  (wave  growth  potential)  and  waves  near  the  shore  (combined  growth 
and  transformation). 

The  wave  exposure  index  developed  by  Keddy  has  been  applied  to  seagrass 
habitats  by  Fonseca  et  al.  (2002).  Although  their  relative  wave  exposure 
index  (REI)  was  a  good  predictor  of  seagrass  spatial  distribution  in  some 
estuaries  (Murphey  and  Fonseca  1995;  Fonseca  and  Bell  1998),  problems 
were  encountered  in  other  areas  (i.e.,  Chiscano  2000).  These  were  likely 
due  to  complex  bottom  bathymetry  such  as  sills  and  sandbars.  As  a  result, 
a  further  refinement  of  the  wave  exposure  index,  the  relative  wave  expo¬ 
sure  index  now  includes  water  depth  and  tidal  fluctuations  (Robbins  et  al. 
2001): 


REI  =  jy,xP,x{idwF,xT)] 

i=l 


(12) 


where  an  inverse  distance  weighting  (IDW)  function  is  used  to  weight  ef¬ 
fective  fetch  (F)  by  depth  from  the  point  of  interest  to  the  shoreline.  Effec¬ 
tive  fetch  is  then  further  weighted  by  the  period  of  tidal  emersion  (T)  de¬ 
rived  from  a  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
tidal  prediction  based  on  tide  curves  generated  from  harmonic  analysis  at 
gauge  locations  (readily  available  on  NOAA’s  Web  site). 

A  spatially  explicit  version  of  the  relative  wave  exposure  index,  the  Wave 
Exposure  Model  (WEMo)  has  been  developed  (Eonseca  and  Bell  1998; 
Kelly  et  al.  2001;  Eonseca  et  al.  2002).  WEMo  calculates  REI  values  at 
specified  locations  as  a  means  of  predicting  the  relative  amount  of  wave 
exposure  at  each  site.  The  model  requires  the  user  to  input  four  variables: 
a  bathymetric  grid  (ARC-based),  a  shoreline  (shapefile),  wind  data  (dura¬ 
tion  and  speed  from  the  eight  major  compass  headings),  and  a  file  that  de¬ 
fines  georeferenced  points  or  sites  of  interest.  Model  output  includes  a  ta¬ 
ble  that  lists  the  REI  value  for  each  point,  a  shapefile  of  the  points,  and  a 
contour  plot  of  interpolated  REI  values.  A  recent  calibration  of  WEMo 
results  using  wave  height  data  collected  in  Chesapeake  Bay  suggest  that 
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WEMo  is  a  good  predictor  of  wave  exposure  in  shallow  seagrass  habitats 
especially  when  25  and  50  percent  exceedance  winds  are  considered.  At 
this  writing,  WEMo  is  being  modified  to  repair  a  dimensionality  problem, 
replace  the  IDW  function  with  published  wave  decay  functions  as  well  as 
adding  user-controlled  selection  of  bottom  friction  values  in  order  to  em¬ 
bed  the  effect  of  seagrasses  on  waves  following  relationships  described  by 
Eonseca  and  Cahalan  (1992).  The  model  will  be  available  for  download 
from  http://www.ccfhr.noaa.gov/  in  late  2006. 

Quantitative  methods 

Wave  height  and  wave  period  are  the  most  desirable  parameters  when 
studying  waves  in  seagrass  habitats.  In  order  to  understand  sedimentary 
processes  in  vegetated  areas  it  is  also  desirable  to  estimate  the  horizontal 
orbital  velocity  at  the  sediment  surface.  These  parameters  can  be  obtained 
through  indirect  calculations  based  on  wind  data  (this  section)  or  direct 
measurements  using  wave  sensors  (discussed  later).  Caution  is  necessary 
when  estimating  wave  parameters  based  on  wind  speed,  fetch  and 
bathymetry  as  other  parameters  such  as  the  presence  of  seagrasses,  mussel 
beds,  and  oyster  reefs,  etc.,  can  alter  the  wave  climate.  Therefore,  these 
calculations  should  be  used  as  estimates  only  and,  whenever  possible, 
should  be  calibrated  against  field  wave  data. 

Estimating  wave  parameters  based  on  wind  and  fetch  data 

To  evaluate  the  success  of  Zostera  marina  transplants  in  the  Dutch  Wad- 
den  Sea  (van  Katwijk  and  Hermus  2000),  significant  wave  height  (average 
of  the  highest  one-third  waves,  Hg)  and  period  (T)  were  estimated  using 
the  Bretscheider  method  (SPM  1984): 


gH. 


u 


2 


0.283  tan 

0.530 

fgd'^ 

0.75  ^ 

tanh 

Ku  / 

0.00565^ 

u 


tanh 


0.530 

fgd^ 

0.75 

(13) 
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=  7.54  tan 

0.833 

'  7 

0.375  ] 

tanh 
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K.U  j 
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0.333 


tanh 


VM  ) 


0.833 

(  gd^ 

1  9 

0.375 

Ku  j 

(14) 


where  g  is  the  acceleration  of  gravity;  u  is  the  wind  velocity  (m/s  at  lo  m 
above  the  water),  d  is  the  average  depth  along  the  fetch  (m),  and  F  is  the 
fetch  (m). 


This  method  takes  into  account  the  bathymetry  of  the  general  area  but 
only  in  the  form  of  an  average  depth  along  the  fetch.  Therefore,  in  areas 
with  sandbars  or  coral/oyster  reefs  this  calculation  overestimates  wave 
parameters  as  the  depth  of  the  sandbar  and/or  reefs  is  likely  to  be  only  a 
small  area  of  the  fetch  being  considered. 


Seagrass  transplants  can  be  compromised  in  areas  of  extensive  sediment 
transport;  therefore,  an  estimate  of  the  horizontal  orbital  velocity  (CT)  at 
the  sediment  surface  may  be  useful  when  selecting  a  potential  restoration 
site.  Van  Katwijk  and  Hermus  (2000)  estimated  this  parameter  for  aZos- 
tera  marina  transplant  site  using  the  formula  proposed  by  Visser  (as  de¬ 
scribed  in  Verhagen  and  van  der  Wegen  1998): 


U  = 


//,  2n:  1 

2  T  sinh(A:J) 


(15) 


where  L  =  -  d/Lo)T,  g  is  the  acceleration  of  gravity,  d  is  the  aver¬ 

age  depth  along  the  fetch,  and  Lo  =  1.56  7^.  Once  again,  caution  is  advised 
as  this  calculation  does  not  take  into  account  local  characteristics  of  the 
site  that  may  alter  the  horizontal  orbital  velocity. 

Measuring  wave  parameters  in  seagrass  habitats 

When  measuring  waves  in  seagrass  habitats,  the  question  to  be  addressed 
needs  to  be  kept  in  mind.  Is  it  necessary  to  determine  the  direction  from 
which  the  waves  are  approaching  the  seagrass  bed  or  will  it  suffice  to  de¬ 
termine  wave  height  and  period?  Wave  direction  data  may  be  relevant 
when  sediment  transport,  pattern  of  seagrass  bed  expansion,  seed  disper¬ 
sal  or  landscape  characteristics  are  being  considered.  Wave  direction 
is  critical  for  estimating  the  magnitude  and  direction  of  wave-driven 
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longshore  currents  and  the  potential  for  seagrass  beds  to  protect  sedi¬ 
ments  and  shorelines  behind  them.  In  contrast,  nondirectional  wave  data 
may  suffice  when  the  questions  posed  address  wave  attenuation,  local 
sediment  resuspension,  etc. 

Quantifying  nondirectionai  waves  in  seagrass  systems 

Visual  estimations.  As  described  in  Morgan  (2000),  a  simple  and  inex¬ 
pensive  way  of  quantifying  waves  is  to  make  a  visual  estimation  of  wave 
height  (vertical  distance  between  troughs  and  crests)  and  period  (time  it 
takes  two  crests  or  troughs  to  pass  a  certain  point)  using  a  vertical  refer¬ 
ence  pole.  This  method  may  be  especially  useful  in  areas  with  oceanic 
waves  (long  periods),  but  in  areas  with  small  and  high-frequency  waves 
such  as  those  often  found  in  estuaries,  this  method  may  be  less  reliable 
(Koch  and  Verduin  2001).  To  aid  in  the  data  gathering  process,  a  video 
camera  focused  on  the  reference  pole  could  be  used  and  later  analyzed  at  a 
slower  speed.  Even  so,  the  data  obtained  using  this  method  are  limited. 

Deployment  of  wave  sensors.  The  use  of  wave  gauges  is  recom¬ 
mended  when  quantifying  wave  parameters  in  seagrass  beds  where  waves 
can  often  be  relatively  small  in  height  and  occur  at  relatively  high  frequen¬ 
cies.  Wave  sensors  are  capable  of  recording  many  times  per  second.  These 
data  are  then  analyzed  to  obtain  the  relevant  wave  parameters. 

Selecting  a  wave  gauge.  There  are  a  series  of  aspects  that  need  to  be 
considered  when  selecting  a  wave  gauge: 

a.  Acoustic  versus  pressure  sensors  -  Many  wave  measuring  instru¬ 
ments  currently  on  the  market  use  acoustic  techniques  (Acoustic 
Doppler  Velocimeter,  acoustic  profilers)  while  others  use  pressure 
transducers.  The  acoustic  techniques  are  less  desirable  as  seagrass 
leaves  interfere  with  the  acoustic  beams  (Koch  and  Verduin  2001). 
Even  so,  acoustic  instruments  can  and  have  been  used  in  seagrass 
beds  (e.g.,  Gacia  et  al.  1999  deployed  at  12-m  depth;  Verduin  and 
Backhaus  2000  worked  at  15-  to  20-m  depth)  but  the  area  where 
measurements  are  being  taken  needs  to  be  cleared  of  seagrass 
leaves.  This  may  lead  to  an  alteration  of  the  local  hydrodynamic 
conditions.  Additionally,  seagrass  debris  and  associated  organisms 
such  as  fishes  can  still  interfere  with  the  acoustic  beam.  Eiltering 
the  data  to  eliminate  these  interferences  is  necessary  when  using 
acoustic  instruments. 
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b.  Water  depth  -  Most  wave  gauges  were  developed  for  oceanic  pur¬ 
poses.  As  a  result,  the  pressure  range  of  these  gauges  is  often  less 
suitable  for  shallow  (<  3  m)  seagrass  habitats.  For  example,  for 
work  in  shallow  water,  the  oceanic  instruments  often  offer  pressure 
transducers  in  the  o  to  30  PSI  range.  In  shallow  seagrass  habitats,  a 
pressure  transducer  in  the  o  to  5  PSI  or  o  to  10  PSI  would  be  more 
appropriate.  Additionally,  some  instruments  stop  recording  when 
exposed  to  air,  therefore,  the  size  of  the  instrument  and  the  de¬ 
ployment  depth  should  also  be  considered.  The  size  of  wave  sensors 
has  decreased  in  recent  years  as  electronic  components  are  becom¬ 
ing  smaller  and  alternative  methods  for  shallow-water  deployment 
are  now  available.  For  example,  in  order  to  be  suitable  for  shallow- 
water  deployment.  Acoustic  Doppler  Velocimeters  and  profilers  can 
now  be  purchased  with  side  looking  probes  such  that  the  instru¬ 
ment  can  be  laid  on  the  sediment  surface  thereby  occupying  a 
smaller  fraction  of  the  water  column. 


c.  Sampling  rate  -  Oceanic  waves  are  longer  than  estuarine  and  shal¬ 
low-water  waves.  Therefore,  a  relatively  low  sampling  rate  (1  to 
2  Hz)  is  suitable  to  describe  oceanic  waves.  In  contrast,  higher  sam¬ 
pling  rates  (2  to  5  Hz)  are  required  to  describe  shallow-water  waves. 
Therefore,  caution  is  needed  when  selecting  an  instrument  to  re¬ 
cord  the  short  waves  often  observed  in  seagrass  habitats.  As  illus¬ 
trated  in  (Figure  8),  the  use  of  an  instrument  that  uses  a  sampling 
interval,  which  is  too  large  to  resolve  the  relatively  short  waves  in 
seagrass  habitats,  may  lead  to  erroneous  data.  The  sampling  rate 
should  be  such  that  at  least  four  samples  are  collected  per  cycle  for 
the  shortest  wave  present  (Denman  1975).  The  highest  frequency 
that  can  be  resolved  at  a  given  sampling  rate  v  is  referred  to  as  the 
Nyquist  frequency  and  is  defined  as: 


Nyquist 


—  V 


(16) 


Note  that  this  is  two  times  higher  than  the  highest  resolvable  wave 
frequency  according  to  Denman  (1975).  For  example,  in  the  Chesa¬ 
peake  Bay,  waves  in  seagrass  habitats  (i-m  depth)  frequently  have 
periods  of  around  2  sec  and  frequencies  of  0.5  Hz.  Therefore,  a 
sample  rate  of  2  Hz  allows  proper  recording  of  these  waves  (f Nyquist 
=  1  Hz  and  f Denman  =  0.5  Hz).  If  significant  wave  energy  exists  in 
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frequencies  higher  than  the  Nyquist  frequency,  it  will  erroneously 
show  up  in  lower  frequencies  in  the  analysis  (aliasing). 


Figure  8.  An  example  of  how  sampling  rate  affects  the  results  of  wave  sampling  protocols. 
The  solid  lines  represent  long/oceanic  (A)  and  short/estuarine  (C  and  E)  waves.  The  x-axis 
represents  time  and  each  vertical  line  represents  one  sample.  When  the  long  wave  (A)  is 
sampled  at  relatively  low  rate  (e.g.,  1  or  2  Hz),  its  characteristics  are  relatively  well 
represented  by  the  data  collected  (B).  When  a  short  wave  (C)  is  sampled  at  the  same  low 
frequency,  the  data  collected  are  not  a  good  representation  (D)  of  the  actual  waves  (C).  In 
contrast,  if  the  same  short  waves  (E)  are  sampled  at  a  higher  rate,  their  characteristics  are 
much  better  represented  by  the  data  collected  (F). 


d.  Wave  attenuation  -  Wave  kinematics  attenuate  with  depth  (veloci¬ 
ties  and  pressures  are  smaller  at  the  bottom  than  near  the  surface). 
Thus,  any  near-bed  measurements  need  to  be  designed  so  that  the 
inherent  noise  in  the  measurements  does  not  dominate  the  wave 
signal  to  be  measured.  For  bottom-mounted  pressure  gauges,  the 
linear  transfer  function  used  to  convert  the  bottom  measurement  to 
the  surface  is  cosh(2JTd/L)2.  For  wavelengths  (or  periods)  where 
this  transfer  function  exceeds  a  value  of  approximately  lOO  to  looo, 
the  measurement  must  be  truncated  and  an  important  portion  of 
the  wave  energy  will  be  lost  if  the  truncation  period  is  not  at  least 
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half  the  peak  wave  period  or  shorter  (Bishop  and  Donelan  1987; 
Smith  2002).  For  example,  a  bottom-mounted  pressure  gauge 
measurement  must  be  in  2  m  water  depth  or  less  to  capture  3-sec 
waves,  while  10-sec  waves  can  be  measured  in  water  depths  up  to 
approximately  20  m.  The  presence  of  strong  current  can  alter  the 
wavelength  and  further  limit  the  depth  of  measurement. 

In  summary,  careful  consideration  is  needed  when  selecting  a  wave  sensor. 
It  may  be  advisable  to  do  a  site  visit  and  to  identify  the  wave  characteris¬ 
tics  of  the  seagrass  habitats  to  be  studied  prior  to  selecting  an  instrument. 
This  could  be  done  in  collaboration  with  physical  oceanographers  or 
coastal  engineers.  Once  general  depth  and  wave  characteristics  for  the 
study  site  are  available.  Table  1  can  be  used  to  guide  the  user  in  the  in¬ 
strument  selection  process. 

Programming  a  wave  sensor.  Each  wave  sensor  comes  with  its  own 
software.  The  user  programs  the  length  of  the  bursts  to  be  recorded  as  well 
as  the  sampling  frequency.  The  sampling  frequency  needed  to  resolve  the 
waves  to  be  studied  has  been  described.  Therefore,  this  section  will  focus 
on  the  sampling  design. 

Waves  are  normally  not  recorded  continuously  over  extensive  periods  of 
time  (days  or  weeks)  due  to  the  volume  of  data  generated.  For  example,  a 
wave  sensor  recording  at  a  5-Hz  rate  will  generate  4,500  data  points  in 
15  min  and  432,000  points  in  a  day.  Therefore,  wave  sensors  are  normally 
programmed  to  turn  on  for  a  determined  period  of  time  (or  number  of 
samples  to  be  recorded)  and  then  “sleep”  until  they  are  asked  to  record 
again.  These  time  intervals  during  which  the  instrument  is  recording  are 
referred  to  as  bursts.  The  length  of  the  bursts  in  seagrass  research  has  var¬ 
ied  from  5  to  13  min  (Gacia  et  al.  1999;  Verduin  and  Backhaus  2000;  Koch 
2002).  A  rule  of  thumb  for  wave  data  collection  is  to  capture  100  to  200 
dominant  wave  periods  to  get  reliable  statistics  (e.g.,  100  waves  with  a  3- 
sec  period  would  be  5  min).  Bursts  are  recorded  at  intervals  pre¬ 
determined  by  the  user  (Figure  9).  The  interval  at  which  each  burst  is  re¬ 
corded  depends  on  the  amount  of  memory  available  as  well  as  how  often 
the  wave  conditions  are  expected  to  change.  During  storm  events,  a  30- 
min  interval  may  be  recommended  while  during  calm  conditions  a  1-  to 
3-hr  interval  will  suffice. 
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Table  1.  Summary  of  wave  sensors  currently  available  on  the  market. 


Sensor 

Type 

Directional 

Minimum 
depth  of 
deployment 

Maximum 

sampling 

rate 

Examples  of  Web 
links 

Pressure 

transducer 

Pressure 

(strain 

gauge) 

No,  unless 

3  sensors 
are  de¬ 
ployed 

0.5  m 

5  Hz 

httD://www.coastal- 

usa.com/ 

Electro¬ 
magnetic 
current 
meter  w/ 
wave  option 

Pressure  + 

electro¬ 

magnetic 

velocity 

Yes  or  no 
(optional) 

0.3  m 

10  Hz 

httD://www.interoce 

ansvs- 

tems.com/s4oDtions 

.htm 

Vector 

velocimeter 

Pressure 
+  acoustic 
velocity 

yes 

0.1  m 

0.5  m 

64  Hz 

httoV/www.nortek- 

as.com/hardware.DhD 

httD://www.sontek.c 

om/oceans.htm 

Current 
meter  and 
profiler 

Pressure 
+  acoustic 
velocity 
profile 

yes 

1  m 

2  Hz 

httD://www.nortek- 

as.com/hardware.DhD 

httD://www.sontek.c 

om/oceans.htm 

Acoustic 
wave  and 
current 
profiler 
(AWAC) 

Acoustic 

surface 

tracking 

+  acoustic 

velocity 

profile 

yes 

2  m 

4  Hz 

httD://www.nortek- 

as.com/hardware.DhD 

httoV/www.rdinstrume 

nts.com/waves.html 

httD://www.sontek.c 

om/oceans.htm 

time 


Figure  9.  Schematic  diagram  of  wave  sensor  sampling  commonly  used  to  record  waves  in 
seagrass  habitats.  The  upper  line  represents  a  routine  in  which  a  burst  (5  to  13  min)  is 
recorded  every  30  min  while  the  lower  line  represents  a  routine  where  a  burst  is  recorded 

every  60  min. 


Deploying  wave  sensors  in  seagrass  habitats.  When  quantifying 
waves  in  seagrass  habitats,  special  attention  needs  to  be  given  to: 
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a.  Location  of  the  instrument  -  the  objective  of  the  wave  sensor  de¬ 
ployment  needs  to  be  kept  in  mind.  If  the  goal  is  to  characterize  the 
waves  that  reach  a  seagrass  bed,  the  wave  sensor  should  be  de¬ 
ployed  seaward  of  the  seagrass  bed  where  waves  are  not  yet  influ¬ 
enced  by  the  vegetation.  In  contrast,  if  the  goal  is  to  characterize  the 
waves  in  a  seagrass  bed,  the  sensor  can  be  placed  directly  in  the  sea¬ 
grass  bed.  In  each  case,  the  distance  from  the  edge  of  the  bed  to  the 
instrument  should  be  recorded  as  well  as  the  shoot  density  and 
canopy  height.  These  data  will  be  useful  in  the  interpretation  of 
data  collected  over  several  seasons  and  for  different  seagrass  spe¬ 
cies.  For  deployment  of  more  than  one  wave  sensor,  see  following 
sections  on  quantifying  direction  waves  and  wave  attenuation. 

b.  Stability  of  the  instrument  -  it  is  essential  that  the  sensor  be  se¬ 
curely  attached  to  a  solid  structure.  Even  the  slightest  motion  or  vi¬ 
bration  can  cause  false  “waves”  in  the  record.  This  can  be  accom¬ 
plished  by  attaching  the  sensor  to  a  tripod  above  the  sediment 
surface  or  pipes  or  poles  pounded  or  jetted  into  the  sediment. 

c.  Interference  with  vegetation  -  if  an  acoustic  instrument  is  used,  the 
seagrasses  need  to  be  cleared  from  the  sampling  volume/path  such 
that,  even  when  the  leaves  bend  over,  they  do  not  interfere  with  the 
acoustic  signal.  Therefore,  the  radius  of  the  area  cleared  should  be 
equal  or  larger  than  the  maximum  seagrass  leaf  length.  This  process 
is  not  necessary  when  using  pressure  transducers  to  quantify  waves. 

Data  processing.  When  purchasing  wave  sensors,  the  manufacturer  of¬ 
ten  provides  software  options  for  wave  analysis.  If  that  is  not  the  case, 
Matlab  ( http://www.mathworks.eom/I  can  be  used  to  process  wave  data.  The 
analysis  required  to  obtain  relevant  wave  parameters  such  as  significant 
wave  height  and  wave  period  is  called  Fast  Fourier  transformation.  This 
analysis  transforms  a  complex  wave  signal  into  a  series  of  sinusoidal  func¬ 
tions  representing  different  wave  frequencies  (Figure  lo).  The  amount  of 
energy  at  each  of  these  frequencies  is  then  plotted  in  the  form  of  wave 
spectra  (Figure  ii).  The  dominant  (peak)  wave  frequency  is  then  used  to 
estimate  the  dominant  wave  period  and  the  integral  of  the  spectrum  (wave 
energy)  is  used  to  calculate  wave  height.  When  waves  propagate  as  swell, 
the  peak  is  generally  well  defined  and  relatively  narrow,  while  when  the 
sea  is  choppy,  the  peak  is  less  well  defined  and  broader  (Figure  ii).  There¬ 
fore,  spectral  width  is  also  an  interesting  parameter  to  obtain.  Addition- 
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ally,  it  is  useful  to  plot  the  wave  spectra  instead  of  just  analyzing  the  nu¬ 
merical  parameters  generated  by  the  program.  When  two  dominant  wave 
frequencies  occur,  double  peaks  are  observed  in  the  spectra.  Generally, 
software  only  analyzes  the  highest  peak  but  the  secondary  peak  may  be 
biologically  relevant.  For  example,  while  wind  waves  in  the  seagrass  habi¬ 
tat  may  generate  one  peak  in  the  spectra,  the  seagrass  leaf  flapping  fre¬ 
quency  may  generate  a  second,  smaller  peak.  The  spectra  can  also  provide 
a  check  on  data  quality,  e.g.,  an  increasing  high-frequency  tail  on  a  spec¬ 
trum  could  indicate  sensor  noise  overwhelming  the  signal  (cutoff  too  high) 
or  unresolved  higher  frequency  waves  (cutoff  too  low). 


Figure  10.  Visual  representation  of  the  wave  breakdown  process  undertaken  in  the  Fast 
Fourier  analysis.  The  original  complex  signal  (bold  line  on  top)  is  decomposed  into  a  series  of 
sinusoidal  functions  representing  different  wave  frequencies  (thinner  lines).  The  amount  of 
energy  in  each  wave/frequency  is  then  plotted  in  the  form  of  a  wave  spectrum  (Figure  11). 

Quantifying  directionai  waves  in  seagrass  systems 

To  measure  directional  waves  requires  a  combination  of  at  least  three  si¬ 
multaneous  measurements  of  pressure,  velocity,  sea-surface  slope,  sea- 
surface  vertical  velocity,  sea-surface  vertical  acceleration,  or  sea-surface 
elevation.  There  are  various  combinations  of  these  measurements  that  can 
provide  the  necessary  information  to  calculate  directional  wave  spectra; 
typical  examples  are  pressure  and  two  velocity  components  at  a  point  loca¬ 
tion  as  with  a  2  axis  current  meter  equipped  with  a  pressure  sensor  (PUV 
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meter),  three  (or  more)  simultaneous  pressure  measurements  that  are 
separated  in  space  by  precisely  known  geometry  as  with  a  pressure  sensor 
array,  or  sea-surface  slope,  vertical  velocity,  and  vertical  acceleration  as 
with  a  wave-rider  buoy. 
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Figure  11.  Nondirectional  wave  spectra  for  two  wave  bursts  recorded  on  the  eastern  side  of 
the  Virginia  portion  of  Chesapeake  Bay  between  Nassawadox  Creek  and  Hungars  Creek.  The 
data  were  coiiected  with  a  5-MHz  Acoustic  Doppier  Veiocimeter  (ADV),  equipped  with  a 
pressure  gauge,  and  mounted  to  a  bottom  ianding  tripod  depioyed  in  3-m  water  depth.  Burst 
number  475  (biue  iine)  shows  a  sharp  peak  at  0.2  Hz  with  a  sharp  deciine  in  energy  at  higher 
frequencies,  representing  a  narrower,  weii-defined  wave  train.  Burst  number  599  (red  iine) 
shows  a  flatter  response  with  a  broad  peak  and  iess  faii-off  of  energy  at  higher  frequencies, 
representing  a  more  chaotic,  choppy  sea  state. 


-  130ct-2005  11:00:11  UTC 

Burst  #  =  475 
Hsig  =  1.0038 
Tp=  5.12 
DTp  =  332 
Dp  =  334 

160ct-2005  01:00:11  UTC 
Burst  #  =  599 
Hsig  =  0.19633 
Tp=  3.5556 
DTp  =  4 
Dp  =  338 


Recently,  Acoustic  Doppler  Current  Profilers  (ADCP)  have  become  popular 
tools  to  measure  directional  waves.  These  instruments  utilize  spatially 
separated  velocity  measurements,  acoustic  beam  surface  tracking,  and 
pressure  data  all  from  a  single  instrument. 

Measuring  directional  waves  in  seagrass  beds  presents  some  challenges 
that  limit  the  options  for  making  these  measurements.  Many  of  the  cur¬ 
rently  popular  techniques  for  measuring  waves  and  currents  involve 
acoustic  instruments,  but  they  do  not  perform  well  in  shallow  seagrass 
beds  due  to  interference  of  the  seagrasses  themselves  with  the  acoustic 
beams  of  the  instruments,  and  also  due  to  the  shallow-water  depth,  often 
less  then  2  m  in  turbid  waters.  PUV  gauges  using  electromagnetic  velocity 
sensors  are  perhaps  one  alternative  for  measuring  directional  waves  in 
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seagrass  beds,  but  Howell  (1998)  described  a  technique  using  Short  Base 
Line  Arrays  (SBLA)  with  three  pressure  sensors  mounted  in  an  equilateral 
triangle  configuration  that  appears  to  be  a  better  option.  This  technique 
offers  many  advantages  such  as  reasonable  cost,  low  maintenance,  and 
well  established  analysis  procedures.  These  SBLAs  have  sensor  spacings 
(Howell  used  1.8  m)  that  are  significantly  shorter  than  a  wavelength,  even 
when  used  in  shallow  water  with  high  frequency  waves  (e.g.,  peak  pe- 
riod=2  sec,  depth  =  1  m,  wavelength  =  6.3  m).  This  feature  allows  calcula¬ 
tions  of  sea-surface  slopes  which  offer  additional  analysis  options.  One  re¬ 
quirement  is  that  the  pressure  sensors  used  on  these  arrays  be  of  high 
quality  and  have  a  resolution  of  0.001  decibars  and  an  absolute  accuracy  of 
0.01  decibars. 

Howell  (1998)  thoroughly  describes  the  computations  and  equations  nec¬ 
essary  to  calculate  directional  wave  parameters  from  SBLA  data  and  offers 
useful  suggestions  on  how  to  reduce  data  storage  requirements  and  pro¬ 
long  deployments.  For  noncommercial  use  there  is  a  publicly  available 
tool,  “DIWASP,  a  directional  wave  spectra  toolbox  for  MATLAB®”  (see 
Figure  12  for  an  example  of  the  output),  that  offers  three  of  the  most 
widely  used  contemporary  analysis  methods  for  directional  waves.  The 
DIWASP  toolbox  and  User’s  Manual  can  be  downloaded  at: 
http://www.cwr.uwa.edu.au/~iohnson/diwasp/diwasp.htnil. 

Quantifying  wave  attenuation  by  seagrasses 

Wave  attenuation  in  kelp  or  seagrass  beds  is  usually  determined  by  meas¬ 
uring  waves  before  they  reach  the  vegetation  and  after  they  have  passed 
through  a  certain  distance  of  vegetated  bottom  (Figure  13).  Therefore, 
when  quantifying  wave  attenuation,  at  least  two  wave  sensors  are  needed. 
If  more  sensors  are  available,  a  series  of  them  can  be  deployed  along  a  line 
perpendicular  to  the  shore  going  through  the  seagrass  bed  and  to  an  area 
offshore  of  the  seagrass  bed  (Figure  13).  This  will  allow  wave  attenuation 
to  be  quantified  for  different  areas  of  the  seagrass  bed:  edge,  center,  etc. 
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Figure  12.  Polar  plot  of  directional  spectrum  for  burst  number  475  (see  Figure  11  for 
additional  information).  Directional  spectrum  was  computed  with  “DIWASP,  a  directional  wave 
spectra  toolbox  for  MATLAB®",  using  the  Iterative  Maximum  Likelihood  Method  (IMLM). 
Radial  distances  indicate  wave  period  in  seconds  and  color  contours  are  power  spectral 
density  which  is  a  measure  of  the  energy  content  coming  from  a  certain  direction  (in  this 
case,  mainly  from  the  NE  and  N  but  also  NW).  For  location  and  methods  see  legend 

of  Figure  11. 

Waves  are  not  only  attenuated  by  the  vegetation  but  also  over  depth.  As  a 
result,  it  is  recommended  that  the  wave  sensors  be  deployed  at  the  same 
depth.  This  can  be  accomplished  by  installing  poles  or  placing  tripods  at 
the  sampling  site  and  attaching  the  wave  gauge  at  the  appropriate  height. 

If  that  is  not  possible,  a  depth  correction  will  be  required.  Huber  (2003) 
accounted  for  wave  shoaling  in  a  Z.  marina  bed  by  applying  the  following 
equation: 


(17) 


where  Ai  and  A2  are  the  wave  amplitudes  (A=H/2.)  at  the  deep  and  shallow 
sites,  respectively;  Cgi  and  Cg2  are  the  group  velocities  at  the  deep  and  shal¬ 
low  sites,  respectively  (Equation  8). 

The  shoreward  edge  of  seagrass  beds  is  usually  defined  by  increasing 
depth.  Therefore,  the  deployment  of  a  wave  sensor  offshore  of  the  seagrass 
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bed  at  the  same  depth  of  the  sensor  deployed  in  the  seagrass  bed  may  be 
difficult.  An  alternative  may  be  to  deploy  one  wave  sensor  in  the  seagrass 
bed  and  one  in  an  adjacent  unvegetated  area  which  has  the  same  depth 
and  fetch  (Figure  13). 


In  order  to  appropriately  interpret  the  results,  the  distance  from  the  wave 
sensor  to  the  edge  of  the  seagrass  bed  should  be  measured  in  different  di¬ 
rections  (all  possible  incident  wave  directions).  If  resources  are  available, 
it  may  also  be  useful  to  determine  seagrass  canopy  height  and  shoot  den¬ 
sity  in  the  area.  This  will  allow  for  interpretation  of  wave  attenuation  as  a 
function  of  these  seagrass  parameters. 


Figure  13.  Alternative  ways  of  deploying  wave  sensors  when  quantifying  wave  attenuation  by 
seagrasses.  The  dark  area  adjacent  to  Bishop’s  Head  Point,  MD,  shows  a  Ruppia  maritima 
bed.  The  left  panel  gives  an  example  of  how  wave  sensors  (red  ovals  and  yellow  stars)  can  be 
deployed  along  a  shore-perpendicular  transect.  Ideally,  all  sensors  will  be  at  the  same  depth. 
The  right  panel  gives  an  example  of  how  wave  sensors  can  be  deployed  at  the  same  depth. 
Red  ovals  represent  essential  wave  gauges  while  yellow  stars  represent  suggested  locations 

for  additional  sensors. 


Wave  sensors  should  be  programmed  to  record  waves  concurrently.  The 
data  can  then  be  plotted  as  wave  height  offshore  versus  wave  height  in  the 
seagrass  bed  (Figure  14).  The  points  usually  falls  below  the  1:1  line  indicat¬ 
ing  that  there  is  wave  attenuation.  The  slope  of  the  regression  line  indi¬ 
cates  the  degree  of  attenuation  of  waves  by  the  seagrass  bed  (Figure  14).  In 
addition  to  plotting  the  data,  it  should  also  be  processed  by  determining 
wave  attenuation  (AT)  as  the  percentage  of  the  original  (offshore)  wave 
height  (Ho)  still  observed  in  the  seagrass  bed  (Hs): 
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Figure  14.  Wave  attenuation  in  a  seagrass  [Ruppia  maritima)  bed  off  Bishop’s  Head,  MD,  in 
June  (top)  and  October  (bottom).  In  June,  shoot  density  was  1,270  shoots  and  plants  were 
reproductive  (occupied  the  entire  water  column)  while  in  October  the  shoot  density  was  1,968 
shoots  m-2  and  plants  were  not  reproductive  occupying  only  a  fraction  of  the  water  column. 
The  black  line  represents  the  1:1  slope,  i.e.,  no  wave  attenuation.  Note  the  higher  wave 

attenuation  in  June  than  in  October. 
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AT  =  ^  =1=100  (18) 

The  ratio  of  the  group  velocities  accounts  for  shoaling  due  to  differences  in 
depth  at  the  observation  points  (if  depths  are  significantly  different, 
changes  in  wave  height  due  to  refraction  should  also  be  considered).  The 
results  are  usually  expressed  as  a  percentage  of  attenuation  over  a  distance 
(e.g.,  percent/m  or  percent/ioo  m).  Wave  attenuation  should  be  deter¬ 
mined  for  each  burst  and  then  plotted  as  a  time  series  of  wave  attenuation 
over  time  (Figure  15). 

The  principles  described  also  apply  to  the  quantification  of  waves  in  wave 
tanks.  Caution  is  required  in  the  wave  generation  process  though.  The 
waves  generated  in  the  lab  should  be  relevant  to  those  observed  in  sea- 
grass  beds.  For  example,  if  a  wave  generator  is  only  able  to  generate  i-sec 
waves  while  local  seagrasses  are  usually  exposed  to  3-  or  4-sec  waves,  the 
attenuation  obtained  in  the  lab  may  not  be  relevant  to  the  field  popula¬ 
tions. 


Figure  15.  Wave  attenuation  in  a  seagrass  (Ruppia  maritima)  bed  off  Bishop’s  Head,  MD,  in 
June  when  piants  were  reproductive  and  occupied  the  entire  water  coiumn.  Note  that  wave 
attenuation  is  a  function  of  water  depth:  at  high  tide  wave  attenuation  is  iow  and  at  iow  tide 

wave  attenuation  is  high. 
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6  Modeling  Waves,  Flow,  and  Sediment 
Transport  in  Seagrass  Systems 

Existing  modeis 

There  have  been  relatively  few  modeling  studies  of  interactions  between 
waves,  flow,  and  sediments  in  seagrass  beds.  Teeter  et  al.  (2001)  provide  a 
comprehensive  review  of  the  physical,  biological,  and  sedimentological 
complexities  involved  in  any  such  effort  and  present  relevant  equations  as 
a  point  of  departure,  but  they  do  not  present  a  complete  model.  They  state 
that  the  primary  limitations  on  developing  a  complete  wave-flow-seagrass- 
sediment  model  are  computational  power  and  information  on  the  fric¬ 
tional  damping  and  bottom  sheltering  effects  of  seagrass  beds.  In  the  face 
of  continuing  improvements  in  computing  power,  frictional  damping  and 
bottom  sheltering  in  seagrass  beds  remain  as  the  key  issues. 

In  fact,  wave/flow  damping  by  aquatic  vegetation  has  been  the  focus  of 
several  previous  modeling  studies.  These  studies  have  focused  on  devel¬ 
opment  of  expressions  and  parameterizations  for  1-  or  2-D  frictional  drag, 
in  terms  of  a  vegetation  Reynolds  number  and/or  vegetation  density.  The 
drag  force  of  the  vegetation  on  the  flow  is  usually  expressed  as: 

F  =  ^pfaU^  (19) 

where  F  is  force  per  volume,  a  is  projected  area  per  volume,  /  is  the  bulk 
drag  coefficient,  and  U  is  the  maximum  wave-induced  velocity  or  flow 
speed.  The  projected  area  a  is  expressed  as  nd  in  Nepf  (iqqqj's  cylinder 
model  for  steady  currents,  where  n  is  the  number  of  shoots  per  unit  bot¬ 
tom  area  and  d  is  a  typical  shoot  diameter.  This  definition  of  a  is  similar 
to  N-b^  in  the  wave  damping  models  by  Kobayashi  et  al.  (1993);  Mendez  et 
al  (1999);  Mendez  and  Losada  (2004);and  Ota  et  al.  (2004)  where  AT  =  n  = 
the  number  of  shoots  per  unit  bottom  area  and  is  defined  as  the  plant 
area  per  unit  height.  Selected  summaries  of  previous  studies  are  pre¬ 
sented  in  the  following  paragraphs. 

Kobayashi  et  al.  (1993)  presented  an  analytical  solution  of  wave  height  de¬ 
cay  through  vegetation  based  on  linear  wave  theory,  a  Reynolds  number 
dependent  drag  parameterization,  and  constant  depth.  They  compared 
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their  model  to  flume  studies  on  artiflcial  kelp  stands  (N  =  iioo  and 
i490/m2,  =  52  mm).  The  flume  studies  consisted  of  60  runs  with  vary¬ 

ing  water  depths  (0.45-0.52  m),  wave  periods  (0.714-2  sec),  and  wave 
heights  (0.036-0.1934  m).  They  used  the  bulk  drag  coefficient  ( / )  to  cali¬ 
brate  the  model  for  60  runs  and  then  correlated  /  with  Reynolds  number 
{R  =  Ud/v,  where  v  is  the  kinematic  viscosity  of  the  water).  They  found 
that  /  decreases  with  increasing  Reynolds  number,  and  the  relationship 
can  be  approximated  by: 


/  =  0.08  + 


r  2200V" 

V  R  J 


(20) 


Equation  20  is  plotted  in  Figure  16a. 


Mendez  et  al.  (1999)  and  Mendez  and  Losada  (2004)  expanded  Kobaya- 
shi’s  solution  by  including  swaying  motion  of  the  seagrass,  wave  breaking, 
and  variable  depth,  and  parameterized  their  model  based  on  careful  flume 
experiments.  They  allowed  for  swaying  motion  of  the  seagrass  by  chang¬ 
ing  the  characteristic  velocity  in  Equation  19  to  the  relative  velocity  be¬ 
tween  plant  and  water.  They  reported  another  empirical  relationship  be¬ 
tween  bulk  drag  and  Reynolds  number: 


/  =  0.4  + 


^4600^ 

V  R  J 


(21) 


Equation  21  is  also  plotted  in  Figure  16a.  Given  the  same  Reynolds  num¬ 
ber,  the  bulk  drag  coefficient  in  Mendez  et  al.  (1999)  is  higher  than  that  in 
Kobayashi  et  al.  (1993)  because  a  lower  velocity  relative  to  the  plant  when 
accounting  for  plant  motion  requires  a  higher  drag  coefficient  to  maintain 
the  same  amount  of  wave  energy  attenuation.  Their  model  fit  to  the  data 
has  a  better  correlation  coefficient  than  Kobayashi  et  al.’s  model,  which 
suggests  that  swaying  motion  of  plant  might  need  to  be  considered  for  op¬ 
timal  drag  estimation. 

Nepf  (1999)  used  a  different  approach  to  explore  the  drag  of  vegetation  on 
steady  currents.  She  ignored  the  flexibility  of  the  vegetation,  mimicking 
the  seagrass  stalks  using  arrays  of  cylinders  (width  d  =  6.4  mm).  Based  on 
observations  for  pairs  of  cylinders  by  Bokaian  and  Geoola  (1984),  she  as¬ 
sumed  that  the  bulk  drag  coefficient  is  a  function  of  vegetation  density  as 
represented  by  the  fractional  volume  occupied  (ad).  Numerical  Simula- 
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tions  were  then  performed  for  both  random  and  staggered  arrays  of  cylin¬ 
ders  with  different  element  spacings  (different  values  of  ad).  She  showed 
that  the  bulk  drag  coefficient  is  relatively  constant  for  ad  up  to  o.oi  and 
declines  steadily  beyond  this  density  (Figure  6  in  Nepf  1999).  An  ap¬ 
proximate  fit  to  her  data  (Equation  25)  is  plotted  in  Figure  16b.  In  the 
density-independent  range  (o.ooi<ad<o.oi),  the  spacing  between  cylin¬ 
ders  is  too  large  for  the  wake  behind  an  upstream  cylinder  to  influence  the 
drag  of  a  downstream  one.  In  the  steady-decline  range  (o.oi<ad<o.i)  the 
drag  coefficient  decreases  due  to  turbulent  wake  interference  that  delays 
the  point  of  separation  on  a  downstream  cylinder  and  subsequently  leads 
to  a  lower  drag  (Kundu  and  Cohen  2002).  It  should  be  noted  that  in  this 
model  /  was  argued  to  be  a  weak  function  of  Reynolds  number. 
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Figure  16.  (a)  Relationship  between  the  bulk  drag  coefficient  ( / )  for  pure  wave  motions  and 
Reynolds  number  (Re).  The  solid  line  (Equation  20)  is  the  relationship  reported  in  Kobayashi 
et  al.  (1993),  whereas  the  dashed  line  (Equation  21)  is  from  Mendez  et  al.  (1999). 

(b)  Approximation  (Equation  25)  of  relationship  between  the  bulk  drag  coefficient  ( )  for 
steady  current  and  the  fractional  volume  occupied  by  vegetation  (aof)  reported  in  Nepf  (1999). 

There  are  two  main  differences  between  drag  models  following  Kobayashi 
et  al.  (1993)  and  the  model  of  Nepf  (1999).  First,  Kobayashi-type  models 
are  for  oscillatory  flow  (waves)  while  Nepf  s  model  is  for  steady  currents. 
Therefore,  the  Reynolds  numbers  are  different  because  the  characteristic 
velocities  (IT)  are  wave  orbital  velocity  and  uniform  current  speed,  respec¬ 
tively.  Second,  in  Kobayashi-type  models  the  bulk  drag  is  a  function  of 
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Reynolds  number  that  reflects  the  nature  of  the  flow  around  a  single  shoot 
of  vegetation.  In  contrast,  the  bulk  drag  in  Nepf  s  model  is  a  function  of 
vegetation  density  that  reflects  a  property  of  the  whole  bed.  The  flume  ex¬ 
periments  presented  in  Kobayashi  et  al.  (1993)  and  Mendez  et  al.  (1999) 
have  an  equivalent  ad  ~  3.  This  value  is  above  the  range  of  ad  (o.ooi~o.i) 
reported  by  Nepf  (1999);  note  that  the  Nepf  values  are  more  representative 
of  natural  seagrass  beds.  The  high  vegetation  density  in  the  Kobayashi  et 
al.  (1993)  and  Mendez  et  al.  (1999)  studies  implies  a  strong  vegetation 
density  dependence  according  to  the  results  of  Nepf  (1999),  although  oscil¬ 
latory  motion  may  have  limited  the  growth  of  turbulent  wakes.  On  the 
other  hand,  the  fact  that  the  Nepf  (1999)  study  was  for  steady  currents  lim¬ 
its  its  direct  applicability  to  wave-seagrass  interactions. 

All  of  these  studies  were  vertically  two-dimensional  (2-D),  measuring  or 
modeling  a  vertical  slice  through  a  seagrass  bed  in  the  direction  of  wave 
propagation  or  flow,  with  flow  prevented  from  diverging  around  the  bed. 
Thus,  while  they  were  all  instructive  and  valuable,  they  could  not  consider 
spatially  varying  seagrass  bed  geometry  (e.g.,  less  than  complete  seagrass 
coverage),  spatially  varying  shorelines  and  bathymetries,  or  combinations 
of  waves  and  currents.  In  addition,  although  there  have  been  several  ob¬ 
servational  studies  that  indicate  enhanced  sediment  deposition  in  seagrass 
beds  (Lopez  and  Garcia  1998;  Gacia  et  al.  1999;  Gacia  and  Duarte  2001; 
Koch  et  al.  2006),  there  have  been  almost  no  modeling  studies  of  sediment 
transport  in  seagrass  beds. 

The  model  described  by  Teeter  et  al.  (2001),  as  implemented  at  least  par¬ 
tially  in  Teeter  (2001),  is  an  exception.  It  is  quite  comprehensive,  includ¬ 
ing  wind  forcing,  wave  forcing,  seagrass-enhanced  drag,  and  sediment 
transport,  but  it  depends  extensively  on  empirical  parameterizations  based 
on  local  observations.  For  example,  Teeter  (2001)  implemented  this 
model  for  Laguna  Madre,  TX,  representing  vegetation  drag  by  a  fixed 
roughness  (/cn  -  0.2  m)  which  was  tuned  to  give  reasonable  agreement 
with  field  observations,  but  is  not  applicable  to  seagrass  beds  in  other  loca¬ 
tions  with  other  combinations  of  waves  and  currents. 

Given  these  shortcomings,  a  new  approach  was  developed  for  modeling 
interactions  between  waves,  currents,  and  sediment  transport  in  seagrass 
systems.  The  new  model  reduces  the  empiricism  of  the  Teeter  et  al.  ap¬ 
proach  by  estimating  wave  and  current  drag  that  depends  on  seagrass 
density  and  height,  based  on  Nepf  (1999).  It  also  considers  three- 
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dimensional  (3-D)  spatial  variability  in  bed  geometry  and  bathymetry,  al¬ 
lows  for  both  wave  and  current  influences,  considers  the  nearshore  cur¬ 
rents  generated  by  wave  breaking,  calculates  total  bottom  shear  stress 
based  on  vector  addition  of  wave  and  current  stresses,  and  estimates  fine 
sediment  resuspension,  deposition,  and  transport  in  and  near  grass  beds. 
In  the  remainder  of  this  section,  the  model  development  is  described  with 
an  emphasis  on  drag  estimation,  validate  it  against  flume  studies  of  flow 
reduction  by  Gambi  et  al.  (1990)  and  against  field  observations  of  wave 
damping  (Koch,  unpublished),  and  present  several  model  scenarios  ex¬ 
ploring  the  effects  of  seagrass  bed  geometries  on  wave  attenuation,  tidal 
current  modification,  and  sediment  trapping. 

Model  development  and  validation 

Numerical  modeling  of  waves  and  currents 

A  quasi-sD,  curvilinear  version  of  the  nearshore  circulation  model 
SHORECIRC  with  wave  driver  REF/DIF  (Kirby  and  Dalrymple  1994;  Shi 
et  al.  2003)  has  been  adapted.  This  modeling  system,  supported  by  the 
U.S.  Nearshore  National  Ocean  Partnership  Program  (NOPP),  aims  to 
predict  current,  wave,  and  wave-driven  current  transformations  in  the 
nearshore  ocean.  SHORECIRC  numerically  solves  the  depth-integrated  2- 
D  horizontal  equations  and  incorporates  a  semianalytical  solution  for  the 
vertical  current  profile  (Svendsen  et  al.  2000).  REF/DIF  accounts  for 
shoaling,  refraction,  energy  dissipation,  and  diffraction  as  waves  propa¬ 
gate  over  variable  bathymetry  and  determines  short-wave  forcing  to  drive 
currents  in  SHORECIRC.  Enhancements  to  the  system  include  estimating 
seagrass  effects  on  drag  and  turbulence,  calculating  the  vector  sum  of  wave 
and  current  bottom  stresses,  and  adding  a  fine  sediment  transport  module. 

For  the  model  presented  here,  the  vegetation  form  drag  expression  of  Nepf 
(1999),  which  was  developed  based  on  laboratory  experiments  with  steady 
flows  through  rigid  grass  mimics,  was  adopted  and  modified.  The  primary 
reason  for  using  this  expression  is  that  it  explicitly  accounts  for  the  effects 
of  seagrass  shoot  density  over  a  realistic  range  of  densities.  The  dominant 
seagrass  species  in  the  field  studies  to  which  model  predictions  was  com¬ 
pared  was  Ruppia  maritima  (leaf  width  -  1.5  mm),  with  a  fractional  vol¬ 
ume  (ad)  that  fluctuated  seasonally  between  about  0.0014  and  0.003.  This 
range  of  ad  is  within  the  density-independent  regime  of  Nepf  (1999),  but  it 
is  three  orders  of  magnitude  smaller  than  the  values  reported  in  Kobayashi 
et  al.  (1993)  and  Mendez  et  al.  (1999)  for  their  laboratory  studies  of  sea- 
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grass  wave  drag.  Steady  flow  drag  data  in  comparable  seagrass  densities 
were  preferable  to  wave  drag  data  from  a  much  higher  seagrass  density, 
especially  since  the  same  basic  drag  formulation  for  both  steady  flow  and 
wave  forcing  in  the  model  was  used. 

Bottom  shear  stress  for  steady  currents  is  written  using  a  standard  quad¬ 
ratic  law: 


T,  =  pC,U^  (22) 

where  p  is  flow  density,  U  is  depth-averaged  flow  velocity,  and  Cd  is  the 
drag  coefficient.  Assuming  that  seagrass  blades  may  be  modeled  as  rigid 
cylinders,  Nepf  (1999)  partitioned  drag  into  skin  friction  due  to  contact 
between  flow  and  sediment  grains  and  form  drag  by  the  seagrass  blades. 
She  expressed  the  drag  coefficient  as: 


1 —  h 

=  (I-  ad)Cg  +  —Cj^ad(—)  (23) 

2  d 

where  a  is  the  projected  plant  area  per  unit  volume,  d  is  shoot  diameter,  h 
is  water  depth,  ad  represents  the  fractional  volume  occupied  by  sea- 
grasses,  Cb  is  a  skin  friction  drag  coefficient  (set  equal  to  0.001  here),  and 
Cjj  is  the  bulk  drag  coefficient  for  seagrass,  which  Nepf  (1999)  determined 
from  experiments.  The  first  term  on  the  right-hand  side  of  Equation  23 
represents  skin  friction,  whereas  the  second  term  represents  form  drag. 
The  calculation  of  a  was  modified  to  allow  seagrasses  to  only  occupy  part 
of  the  water  column,  so  a  =  nld/h,  where  n  is  the  number  of  seagrass 
shoots  per  unit  area  and  /  is  the  shoot  height.  Substituting  into  Equa¬ 
tion  23,  the  drag  coefficient  for  current  becomes: 

nld^  1  — 

Q=(l - —)Cs+-Cijnld  (24) 

n  2 

Thus,  it  is  a  function  of  shoot  height,  shoot  density,  shoot  diameter,  and 
water  depth.  In  Nepf  (i999)’s  model,  the  bulk  drag  coefficient  is  a  func¬ 
tion  of  fractional  volume  (ad).  We  approximate  the  curve  in  Eigure  16b  of 
Nepf  (1999)  as: 


—  1.17,  10’^  <  ad  <  10'" 

C  ~ 

""  -0.255  In(ad),  10‘"  <  ad  <  10"' 


(25) 


ERDCTR-06-15 


43 


In  Nepf  s  experiments  the  velocity  was  measured  at  7.5  cm  above  the  bot¬ 
tom,  whereas  the  reference  height  in  SHORECIRC  was  set  at  1  m.  Hence, 
the  drag  coefficient  for  7.5  cm  height  needs  to  be  converted  to  that  for  1  m. 
Assuming  that  the  near-bottom  velocity  profile  is  logarithmic,  we  calculate 
the  bottom  roughness  coefficient  Zo  using  the  well-known  relationship  be¬ 
tween  the  drag  coefficient  and  the  bottom  roughness  parameter  in  a  rough 
turbulent  boundary  layer  (Equation  26).  Cd  is  obtained  from  Equation  24 
and  reference  height  z  equal  to  7.5  cm;  k  is  the  von  Karman  constant, 
equal  to  0.4. 


-\2 


C 


d 


In 


V^oy 


(26) 


Einally,  using  this  Zq  and  a  new  reference  height  of  1  m,  the  drag  coeffi¬ 
cient  for  SHORECIRC  is  obtained. 

Eield  observations  were  used  to  determine  the  seagrass  bulk  drag  for 
waves  in  REE/DIE.  Bottom  shear  stress  due  to  pure  wave  action  ( is 

expressed  in  terms  of  the  wave  friction  factor  (f^: 

(27) 


where 


/ 


1- 


nld 


\fB+\fhld 


(28) 


where  Ub  is  wave  orbital  velocity  near  bottom, /b  is  the  wave  skin  friction 
factor,  and  /  is  a  bulk  drag  representing  the  effects  of  seagrasses  on 
waves,  /b  was  calculated  using  a  bottom  roughness  equivalent  to  the  value 
of  Cb  =  0.001  used  in  SHORECIRC,  following  procedures  in  U.S.  Army 
Corps  of  Engineers  (2002),  while  /  was  determined  using  field  observa¬ 
tions  and  assuming  the  functional  form  of  Equation  25  with  an  adjustable 
multiplicative  coefficient. 
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Once  the  drag  coefficient  and  wave  friction  factor  are  estimated  through 
Equations  24  and  28,  current  and  wave  fields  are  calculated  by 
SHORECIRC  and  REE/DIE.  With  this  updated  current  and  wave  field 
(wave  height  and  period)  and  with  known  bottom  sediment  grain  size,  skin 
friction  shear  stress  due  to  pure  current  (  and  wave  motions  (  can 
be  obtained  using  the  techniques  in  U.S.  Army  Corps  of  Engineers  (2002). 
Then  we  apply  vector  summation  of  the  two  skin  friction  components  to 
calculate  maximum  skin  friction  shear  stress  (  as  follows: 

(29) 

where  is  the  angle  between  current  and  wave  propagation  and  can  be 
calculated  from  SHORECIRC  and  REE/DIE.  Because  we  are  interested  in 
the  maximum  potential  for  sediment  movement,  the  absolute  value  of 
is  used  in  Equation  29  to  guarantee  maximum  vector  summation 
regardless  of  the  direction  of  wave  orbital  motion.  This  vector  summation 
ignores  enhanced  turbulence  due  to  nonlinear  wave-current  interactions 
in  the  bottom  boundary  layer  (Grant  and  Madsen  1979).  However,  given 
the  high  uncertainty  of  seagrass  drag  estimation  and  turbulence  structure 
in  seagrass  beds.  Equation  29  is  a  reasonable  first  order  approximation  for 
combined  wave-current  bottom  stress. 

Sediment  transport  modeling 

A  suspended  sediment  transport  module  has  been  developed  and  incorpo¬ 
rated  based  on  North  et  al.  (2004).  The  module  accounts  for  erosion  and 
deposition  with  a  simple  parameterization  of  consolidation  for  single¬ 
grain-size  cohesive  sediments.  We  solve  for  changes  in  bottom  sediment 
per  unit  area  (B  in  kg/m^)  over  time  t  at  each  grid  point  using: 

dB 

—  =  D-E-]B  (30) 

dt 

where  D  and  E  are  the  deposition  and  erosion  rate  (kg/m/sec),  respec¬ 
tively,  and  7  is  a  first-order  consolidation  rate  (see  k  set  equal  to  zero 
here).  The  formulation  states  that  the  amount  of  erodible  sediment  per 
unit  area  increases  by  deposition  but  decreases  by  erosion  and  consolida¬ 
tion.  The  deposition  rate  is  calculated  as: 


D  =  W,C 


(31) 
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where  the  settling  velocity  (tVs)  is  equal  to  0.03  cm/sec  (a  typical  value  for 
fine  suspended  sediment  in  Chesapeake  Bay,  Sanford  et  al.  2001)  and  C  is 
depth-averaged  suspended  sediment  concentration  (kg/m3).  The  erosion 
rate  may  be  expressed  as: 


E  =  +  2^)  (32) 

E  dt 

where  r^is  critical  shear  stress  for  erosion  (e.g.,  approximately  0.15  Pa  for 
fine  sand),  M  is  an  empirical  constant  (5x10-5  kg/m/sec  here),  and  H  is  the 
Heaviside  step  function  {H=i  when  its  argument  is  >  o  and  H=o  when  its 
argument  is  <  o).  The  first  step  function  in  Equation  28  represents  the  ini¬ 
tiation  of  sediment  motion  when  the  maximum  bottom  shear  stress  ex¬ 
ceeds  the  critical  value,  while  the  second  step  function  is  to  prevent  over¬ 
erosion  and  negative  values  of  B. 


Given  the  erosion  and  deposition  rates  in  each  model  cell,  a  third-order 
accurate  numerical  scheme  QUICKEST  (Leonard  1979)  is  used  to  solve  the 
depth-averaged  transport  equation  for  suspended  sediments  (Clarke  and 
Elliot  1998): 


d(HC)  ^  d(HUC)  ^  d(HVC) 
dt  dx  dy 


ox  dx  oy  dy 


(33) 


where  H  is  water  depth,  U  and  V  are  depth-averaged  velocity  components, 
and  Kx  and  Ky  are  diffusion  coefficients.  The  QUICKEST  scheme  has  been 
shown  to  avoid  overshoot  problems  near  strong  gradients  in  concentra¬ 
tion.  This  feature  is  particularly  important  because  the  presence  of  sea- 
grass  could  lead  to  abrupt  changes  in  bottom  shear  stress,  which  may  in 
turn  cause  strong  gradients  in  suspended  sediment  concentrations. 

Model  validation 

Eirst,  the  model  to  test  the  effects  of  seagrasses  under  current-only  condi¬ 
tions  (SHORECIRC  only)  was  set  up.  Because  we  are  particularly  inter¬ 
ested  in  modeling  the  effects  of  seagrass  beds  that  cover  only  part  of  the 
model  domain,  such  that  water  may  flow  around  the  bed  rather  than  being 
forced  through  it  or  over  it,  the  data  of  Gambi  et  al.  (1990)  was  used  for 
comparison.  Gambi  et  al.  studied  flow  speed  reduction  by  Zostera  marina 
L.  (eelgrass)  in  a  seawater  flume,  with  the  seagrass  bed  occupying  only 
20  percent  of  the  width  of  the  flume.  SHORECIRC  was  configured  to 
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mimic  the  relative  dimensions  of  their  flume  experiments.  The  actual 
model  domain  was  considerably  larger  because  of  computational  con¬ 
straints,  but  the  ratios  of  the  domain  length,  domain  width,  and  the  size  of 
the  eelgrass  bed  were  scaled  exactly.  Flow  was  driven  using  an  upstream 
flux  boundary  condition,  with  no  flow  through  the  domain  sidewalls,  to 
generate  the  same  free-stream  velocities  as  Gambi  et  al.  The  eelgrass  pa¬ 
rameters  they  reported  were  used  to  calculate  the  drag  coefficient  for 
SHORECIRC  based  on  Equations  24  and  25.  The  volume  flux  reduction 
was  computed  within  the  eelgrass  bed  from  just  upstream  of  the  bed  to  the 
end  of  bed,  the  same  relative  locations  as  in  their  flume  measurements. 

The  volume  flux  reduction  is  defined  as: 

(l-j[/dz/jt/_,„,dz)xlOO  (34) 

where  Ucontroi  is  the  upstream  velocity.  Combinations  of  two  shoot  densi¬ 
ties  (600  and  1,200  shoots/m2)  and  two  free-stream  velocities  (10  and 
20  cm/sec)  were  chosen.  Comparisons  are  shown  in  Eigure  17. 


Figure  17.  Comparison  of  flow  speed  reduction  by  a  seagrass  bed  between  modei  predictions 
(4  curves)  and  flume  data  reported  by  Gambi  et  ai.  (1990)  (4  discrete  points).  The  modei  is 
set  up  using  the  same  reiative  dimension  as  Gambi  et  ai.’s  flume,  and  the  same  seagrass 
parameters  such  as  shoot  density,  shoot  height,  and  ieaf  width  are  appiied.  The  percent  of 
voiume  flux  reduction  is  used  as  an  indicator  of  flow  speed  reduction  and  is  piotted  against 
the  distance  downstream  from  the  free-stream  veiocity  measurements.  The  ieading  edge  of 
the  seagrass  bed  is  indicated  by  the  verticai  iine  at  320  m.  The  open  and  solid  triangies 
represent  shoot  densities  of  l,200/m2  with  20  and  10  cm/sec  free-stream  veiocities, 
respectiveiy,  whereas  the  open  and  soiid  circies  are  GOO/m^  with  20  and  10  cm/sec. 
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The  model-predicted  values  for  the  four  different  combinations  agree  rea¬ 
sonably  well  with  Gambi  et  al.’s  results,  without  any  parameter  tuning.  As 
expected,  the  eelgrass  bed  with  higher  shoot  density  results  in  higher  vol¬ 
ume  flux  reductions.  The  model-predicted  volume  flux  reduction  increases 
rapidly  behind  the  leading  edge  of  eelgrass  bed  and  levels  off  approxi¬ 
mately  halfway  into  the  bed. 

Next,  the  model  was  set  up  to  test  the  effects  of  seagrass  on  wave  attenua¬ 
tion  (RED/DIF  only).  As  mentioned  already,  the  formulation  (Equa¬ 
tion  28)  for  estimating  the  wave  friction  factor  has  a  similar  functional 
form  as  that  for  steady  currents.  However,  equivalent  data  to  that  of  Nepf_ 
(1999)  on  the  relationship  between  seagrass  density  and  wave  form  drag  / 
is  not  available.  Therefore,  field  observations  were  used  to  determine  the 
magnitude  of  the  wave  form  drag.  The  field  observations  were  carried  out 
in  Duck  Point  Cove,  near  Bishop’s  Head  Point,  MD,  in  mesohaline  Chesa¬ 
peake  Bay  (Newell  and  Koch  2004).  Time  series  of  wave  height  and  sea¬ 
grass  parameters  were  measured  in  different  months  at  two  adjacent  sites 
parallel  to  the  shoreline,  one  vegetated  with  R.  maritima  and  the  other 
unvegetated.  The  size  of  R.  maritima  bed  was  about  600  m  in  the  along¬ 
shore  direction  and  200  m  in  the  cross-shore  direction,  and  a  pressure 
sensor  was  located  at  the  center  in  average  water  depth  of  1  m.  Assuming 
the  same  incident  wave  climates  at  the  two  sites,  wave  height  measure¬ 
ments  can  be  plotted  at  the  unvegetated  site  against  the  vegetated  sjte  to 
evaluate  wave  attenuation  by  the  R.  maritima  bed.  Assuming  that  /  is  a 
function  of  fractional  volume  (ad)  and  has  similar  functional  form  to  that 
for  steady  current  (Equation  25),  the  height  of  the  /  curve  was  changed  to 
obtain  the  observed  wave  attenuations  in  October.  The  calibrated  /  is 
written  as: 


-  0.253,  10'' <  ad  <  10"" 

,  (35) 

-  0.055  In(ad),  10'' <  ad  <  10"* 

Applying  observed  seagrass  parameters  for  May  and  June,  the  correspond¬ 
ing  wave  friction  factors  were  calculated  and  the  model-predicted  wave 
height  with  observations  were  compared.  Figures  18  (a),  (b),  and  (c)  show 
the  comparisons  in  May,  June,  and  October,  respectively. 

Wave  attenuation  by  the  seagrass  bed  peaks  in  June  when  the  shoot  height 
occupied  the  whole  water  column  as  a  result  of  the  development  of  repro¬ 
ductive  shoots.  The  slope  of  the  linear  fit  between  the  unvegetated  and 
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vegetated  wave  heights  in  June  was  0.68  (a  slope  of  1  means  no  wave  at¬ 
tenuation).  The  model  qualitatively  captures  this  trend  (June  >  May  and 
Oct)  although  the  model  tends  to  slightly  underestimate  wave  attenuation 
in  June. 


(a)  (b) 


(0 


wave  height  at  unvegetated  site  (m) 

Figure  18.  Comparison  between  model-predicted  wave  attenuation  by  a  seagrass  bed  and 
field  observations.  Each  panel  shows  field  data  (open  circle),  linear  fit  of  field  data  (solid  line 
and  regression  function),  1:1  slope  (dashed  line),  and  model  prediction  (-x-)  for  a  month.  The 
shoot  density  of  Ruppia  maritima  bed  is  about  1536, 1270,  and  1968/m2  for  May  (a),  June 
(b),  and  October  (c),  respectively,  and  the  corresponding  canopy  height  is  0.4, 1.0,  and 
0.4  m.  The  model  Is  calibrated  by  use  of  October  data  (see  text).  The  1:1  line  represents  no 
wave  attenuation.  From  the  linear  regression  of  field  data,  wave  attenuation  Is  higher  in  June 
than  in  May  and  October  due  to  the  presence  of  long  reproductive  shoots. 

This  approach  has  the  advantage  that  a  wider  range  of  vegetation  density 
is  covered  with  one  empirical  parameter  ( /).  This  is  particularly  useful  for 
simulating  seasonally  or  geographically  varying  seagrass  populations.  The 
underestimation  of  wave  attenuation  in  June  may  be  due  to  a  different  re¬ 
sponse  to  oscillatory  forcing,  the  flexibility  of  natural  seagrass  blades,  or  a 
Reynolds  number  dependence  not  accounted  for.  However,  for  the  pre¬ 
sent  purpose  the  fact  that  a  qualitative  reproduction  of  changing  wave  drag 
was  achieved  due  to  seasonal  changes  in  seagrass  morphology  is  consid- 
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ered  to  be  sufficient.  Further  study  is  needed  to  understand  seagrass  bulk 
drag  accounting  for  all  of  these  effects  as  well  as  a  wide  range  of  realistic 
vegetation  densities. 

Model  scenarios 

The  model  domain  is  set  at  720  m  in  the  shore-normal  direction  and 
5,400  m  in  the  shore-parallel  direction  with  a  10-m  x  30-m  grid  resolu¬ 
tion,  respectively  (Figure  19).  Two  bathymetries  are  set  up:  a  flat  bottom 
with  i-m  depth  and  a  sloping  bottom  with  2.5-m  offshore  and  0.05  m  at 
the  shoreline.  When  present,  tidal  currents  are  assumed  to  be  primarily  in 
shore-parallel  direction  with  a  maximum  magnitude  of  about  20  cm/sec. 
Tidal  currents  are  simulated  by  imposing  flux  boundary  conditions 
through  the  upstream  and  downstream  boundaries  of  the  domain  at  semi¬ 
diurnal  frequency.  A  4-s  sinusoidal  wave  enters  the  domain  from  the  off¬ 
shore  boundary  with  wave  heights  varying  between  0.1  and  0.4  m,  at  an 
incident  angle  of  either  o  deg  (Scenarios  1-3)  or  10  deg  (Scenario  4)  coun¬ 
terclockwise  from  the  shore-normal  direction.  The  domain  of  the  sediment 
module  is  smaller  than  the  entire  SHORECIRC/REFDIF  domain  to  avoid 
anomalous  physical  forcing  near  the  boundaries,  and  a  looping  boundary 
condition  is  applied  in  the  shore-  parallel  direction  so  that  the  sediment 
flux  leaving  one  end  of  the  domain  equals  the  flux  entering  the  other  end 
of  the  domain.  Bottom  sediments  are  initialized  with  B  =  3  kg/m^  uni¬ 
formly  distributed  throughout  the  domain.  In  addition  to  the  scenarios 
reported  here,  the  sediment  transport  module  was  verified  to  conserve 
mass  when  suspended  sediments  and  bottom  sediments  were  totaled. 


along-shore  distance  (m) 
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Figure  19.  Model  domain  (top  view)  and  two  bathymetries  (side  view)  for  model  scenarios 
described  in  Table  2.  The  domain  is  720  m  in  cross-shore  direction  and  5,400  m  in  along¬ 
shore  direction  with  a  10-m  x  30-m  rectangular  grid.  The  bathymetry  with  constant  depth 
(1-m)  is  used  in  Scenarios  1  to  3,  and  the  sloping  bottom  with  2.5-m  offshore  depth  and 
0.05-m  onshore  depth  is  used  in  Scenario  4. 
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Table  2.  Model  scenarios. 


Scenario 

Bathymetry 

Cross¬ 
shore  Bed 
Width  (m) 

Along¬ 
shore  Bed 
Length 
(m) 

Position 

Physical 

Forcings 

Output 

Quantities 

1 

Flat(l  m) 

0  to  700 

(50m 

interval) 

Full  along¬ 
shore 
domain 

Offshore 
edge  of  bed 
fixed  at  off¬ 
shore 
boundary 

0.1  to  0.4  m 
waves,  4 
sec,  inci¬ 
dent  angle 
=  0  (shore- 
normal) 

%  of  wave 
energy 
flux  re¬ 
duction 

2 

Flat(l  m) 

100 

300  to 
full  along¬ 
shore 
domain 
(300m 
interval) 

Offshore 
edge  of  bed 
fixed  at  off¬ 
shore 
boundary; 
Center  of 
bed  fixed  at 
the  center 
of  along¬ 
shore 
domain 

Same  as 
above 

Same  as 
above 

3 

Flat(l  m) 

100 

Full  along¬ 
shore 
domain 

Center  of 
bed  located 
50  to  650 
m  from  off¬ 
shore 
boundary 
(50-m 
interval) 

Same  as 
above 

%  of  total 
bottom 
stress 
reduction 

4 

Sloping 
{2.5m  off¬ 
shore, 

0.05m  on¬ 
shore) 

0 

0 

0.1  m,  4 
sec  waves, 
angle  =  10; 
alongshore 
tides 

Wave 

height 

and  skin 

friction 

shear 

stress 

200 

1800 

Center  of 
bed  located 
550  m  from 
offshore 
boundary 

Same  as 
above 

Same  as 
above 

Full  cross¬ 
shore 
domain 

1800 

Center  of 
bed  located 
360  m  from 
offshore 
boundary 

Same  as 
above 

Same  as 
above 
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The  seagrass  parameters  observed  in  June  for  R.  maritima  are  applied 
(density  is  i,27o/m2;  shoot  height  is  i  m).  The  circulation,  wave,  and  sedi¬ 
ment  modules  are  turned  on  in  all  scenarios,  and  we  look  at  several  output 
quantities.  For  the  first  two  scenarios,  the  output  quantity  is  the  percent¬ 
age  of  wave  energy  flux  reduction.  Wave  energy  flux  (F)  is  the  rate  at  which 
wave  energy  is  transported  in  the  horizontal  direction  and  can  be  ex¬ 
pressed  as: 


F  =  EC^  =  (^pgH^)C^  (36) 

where  E  is  the  wave  energy  density,  is  group  velocity,  p  is  water  den¬ 
sity,  g  is  the  gravitational  constant,  and  H  is  wave  height.  In  these  two  sce¬ 
narios  the  geometry  of  the  seagrass  bed  is  changed  and  the  ratio  of  F  is 
calculated  with  and  without  seagrasses,  averaged  over  the  entire  shoreline. 
The  percentage  of  wave  energy  flux  reduction  is  then  (i-F  with/ F without)  X 
100.  For  the  third  scenario,  the  output  quantity  is  the  percentage  of  total 
bottom  stress  reduction.  Because  we  are  interested  in  the  influences  of  the 
seagrass  bed  on  the  total  force  acting  on  the  bottom  sediments  in  the  do¬ 
main,  the  total  bottom  stress  is  defined  as  the  skin  friction  shear  stress  in¬ 
tegrated  over  the  whole  domain.  The  ratio  of  total  bottom  stress  with  and 
without  the  seagrass  bed  is  used  to  calculate  the  percentage  reduction.  In 
the  third  scenario,  the  percent  bottom  stress  reduction  is  compared  as  the 
midpoint  of  the  bed  is  moved  from  an  inshore  position  toward  the  offshore 
boundary,  with  bed  width  and  length  fixed.  The  last  scenario  (sloping  bot¬ 
tom)  is  more  realistic  than  the  constant  depth  scenarios,  and  the  model 
outputs  are  examined  in  more  detail.  Changes  are  examined  in  wave 
height,  skin  friction  shear  stress,  suspended  sediment,  and  bottom  sedi¬ 
ment  over  both  space  and  time  through  two  tidal  cycles.  Figure  20  shows 
changes  in  wave  energy  flux  reduction  when  the  cross-shore  bed  width  is 
varied  but  the  bed  occupies  the  entire  domain  in  the  along-shore  direction 
(Scenario  1). 

The  results  are  presented  with  respect  to  only  the  cross-shore  direction, 
since  there  is  no  alongshore  variation.  Wave  energy  flux  reduction  in¬ 
creases  with  cross-shore  coverage  but  levels  off  as  complete  coverage  is 
approached.  The  increase  in  energy  flux  reduction  is  obviously  due  to  the 
increase  in  seagrass  wave  drag  as  the  bed  becomes  wider.  The  energy  flux 
reduction  levels  off  at  high  percent  coverage  simply  because  not  much 
wave  energy  is  left  to  dissipate,  so  the  rate  of  change  decreases.  Expressed 
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mathematically  using  linear  wave  theory,  the  variation  of  wave  energy  flux 
can  be  written  as: 


dF  of  Q.5HC0. 

—  =  - y 

dx  ^71  kh 


(37) 


where  e  is  energy  dissipation  rate, /is  wave  friction  factor,  co  is  angular 
frequency,  k  is  wave  number,  h  is  water  depth,  and  x  is  distance  across- 
shore.  Percent  energy  flux  reduction  increases  with  increasing  wave 
height  because  a  larger  wave  height  exerts  a  higher  stress  on  the  bottom, 
proportional  to  the  wave  orbital  velocity  squared.  The  wave  energy  dissi¬ 
pation  rate  is  proportional  to  the  product  of  stress  and  wave  orbital  veloc¬ 
ity  for  rough  turbulent  flow  (Dean  and  Dalrymple  1991),  so  wave  energy 
dissipation  is  proportional  to  orbital  velocity  (thus  wave  height)  cubed, 
while  wave  energy  flux  is  only  proportional  to  wave  height  squared.  Thus, 
wave  energy  dissipation  is  proportionately  more  effective  for  higher  waves. 


Figure  20.  Changes  in  wave  energy  flux  reduction  (indicator  of  wave  attenuation)  at  the 
shoreline  when  the  cross-shore  coverage  of  the  seagrass  bed  increases  (see  Scenario  1  in 
Table  2  for  details).  Four  different  incident  wave  heights  (0.1  to  0.4  m)  are  applied  and 
represented  by  different  symbols  (star,  circle,  diamond,  and  cross).  Incident  waves  are  shore- 

normal  direction. 
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Figure  21  shows  changes  in  wave  energy  flux  reduction  for  Scenario  2,  in 
which  the  alongshore  coverage  of  the  bed  is  changed  while  keeping  the 
cross-shore  coverage  fixed.  As  expected,  wave  energy  flux  reduction  is 
linearly  proportional  to  the  alongshore  seagrass  coverage.  Again,  percent 
energy  flux  reduction  increases  with  incoming  wave  height. 

Figure  22  presents  the  percent  reduction  in  total  force  acting  on  bottom 
sediments  as  the  position  of  the  bed  is  moved  from  inshore  towards  the 
offshore  boundary  (Scenario  3),  with  a  fixed  bed  width  of  100  m  and 
length  covering  the  whole  domain  in  the  alongshore  direction. 

It  can  be  seen  that  total  bottom  stress  reduction  increases  approximately 
linearly  with  the  offshore  distance  of  the  bed.  It  makes  sense  that  the  total 
force  acting  on  bottom  sediments  is  reduced  by  moving  the  bed  offshore 
because  the  affected  area  between  the  bed  and  shoreline  increases  linearly 
with  the  distance  of  the  bed  offshore.  The  smaller  waves  that  emerge  from 
the  seagrass  bed  act  over  this  entire  area.  Again,  larger  waves  result  in  pro¬ 
portionately  higher  bottom  stress  reduction. 


%  of  coverage  by  a  100  m-wide  seagrass  bed 
in  along-shore  direction 


Figure  21.  Changes  in  wave  energy  flux  reduction  at  the  shoreline  when  the  alongshore 
coverage  of  the  seagrass  bed  increases  with  fixed  cross-shore  coverage  (Scenario  2).  Four 
different  incident  wave  heights  (0.1  to  0.4  m)  are  applied  and  represented  by  different 
symbols.  Incident  waves  are  shore-normal. 
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For  the  fourth  scenario,  the  more  realistic  sloping  bottom  case,  compari¬ 
sons  are  made  of  model  runs  with  no  seagrass,  a  seagrass  bed  200  m  wide 
and  1,800  m  long,  and  a  seagrass  bed  that  covers  the  entire  width  of  the 
domain  and  is  1,800  m  long.  Figure  23  shows  across-shore  transects  of 
wave  height  and  skin  friction  shear  stress  across  the  center  of  the  seagrass 
bed  at  slack  tide.  In  Figure  23a,  wave  shoaling  and  then  breaking  as  waves 
propagate  shoreward  can  be  seen  without  the  seagrass  bed.  This  wave 
height  evolution  corresponds  to  the  increase  and  quick  drop  of  skin  fric¬ 
tion  shear  stress  shown  in  Figure  23b.  In  both  cases  with  seagrass  beds, 
wave  height  and  skin  friction  shear  stress  within  and  behind  the  bed  are 
greatly  reduced.  The  breaking  zone  and  the  peak  of  skin  friction  shear 
stress  for  the  case  with  a  200-m-wide  bed  are  moved  shoreward.  When 
the  cross-shore  domain  is  fully  occupied  by  the  seagrass  bed,  the  breaking 
zone  disappears.  Differences  between  all  three  cases  in  deeper  water  near 
the  offshore  boundary  are  relatively  small.  The  reason  is  that  short  period 
wave  orbital  velocity  decays  with  depth,  making  bottom  friction  less  effec¬ 
tive  to  dissipate  wave  energy  in  deeper  areas.  Thus,  interactions  between 
seagrass  beds  and  waves  in  deeper  water  depend  on  wave  period;  longer 
period  waves  interact  more  effectively  with  seagrass  beds  in  deeper  water. 


Figure  22.  Changes  in  total  bottom  stress  reduction  (indicator  of  total  force  acting  on  bottom 
sediments)  when  a  seagrass  bed  with  fixed  alongshore  and  cross-shore  coverage  (100  m 
wide)  is  moved  offshore  (Scenario  3).  Four  different  incident  wave  heights  (0.1  to  0.4  m)  are 
applied  and  represented  by  different  symbols.  Incident  waves  are  shore-normal. 
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(a) 


Figure  23.  (a)  wave  height  evoiution  when  waves  propagate  aiong  the  centrai  transect  from 
the  offshore  boundary  to  shoreiine  under  three  seagrass  bed  configurations:  no  seagrass 
(star),  fuli-coverage  (circie),  and  200-m-wide  bed  (diamond).  The  iocation  of  200-m-wide  bed 
is  indicated  by  arrows  and  is  iiiustrated  in  Figure  24  (a)  by  the  rectanguiar  box.  See  Scenario 
4  in  Tabie  2  for  detaiis.  (b)  Changes  in  skin  friction  shear  stress  corresponding  to  wave  height 
evoiution  with  cross-shore  distance  in  (a).  The  soiid  iine  represents  the  criticai  shear  stress 
above  which  fine  sands  of  0.2  mm  diameter  start  to  move. 


Reduced  skin  friction  has  important  implications  from  the  standpoint  of 
sediment  transport.  To  demonstrate  this,  a  line  is  placed  in  Figure  23b  to 
indicate  the  critical  shear  stress  (about  0.15  Pa;  U.S.  Army  Corps  of  Engi¬ 
neers  2002)  for  fine  sands  (0.2  mm).  Sediments  start  to  move  when  shear 
stress  exceeds  a  critical  value  (Equation  32).  As  shown  in  Eigure  23b,  the 
distances  over  which  the  critical  stress  is  exceeded  are  about  the  same  with 
or  without  the  seagrass  beds.  However,  the  erosion  rate  is  proportional  to 
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the  distance  between  the  lines  of  wave-induced  skin  friction  and  critical 
shear  stress  (Equation  32).  Erosion  rate  is  thus  greatly  reduced  within  and 
behind  the  seagrass  beds.  This  implies  that,  without  advection  of  sediment 
from  external  sources,  suspended  sediment  concentrations  within  and  be¬ 
hind  the  beds  may  be  lower  than  those  with  no  seagrass  bed.  Although 
greatly  simplified,  these  model  results  illustrate  that  realistic  bed  geome¬ 
tries  can  have  profound  effects  on  waves  and  can  subsequently  influence 
sediment  dynamics. 

To  further  examine  the  effects  of  seagrass  beds  on  sediment  dynamics,  the 
time  series  is  compared  of  six  variables  associated  with  sediments  between 
the  200-m-wide  bed  case  (Eigure  24a)  and  the  no  seagrass  case  (Eig- 
ure  24b)  over  two  full  tidal  cycles.  The  variables  are  bottom  sediments, 
suspended  sediment  concentration,  skin  friction  shear  stress,  current 
magnitude,  erosion  rate,  and  deposition  rate.  In  Eigure  24,  each  panel  con¬ 
tains  three  lines  that  represent  the  averaged  values  of  each  variable  off¬ 
shore  of  the  bed,  within  the  bed  (or  where  the  bed  would  be),  and  between 
the  bed  and  the  shoreline.  As  can  be  seen  in  Eigure  24,  current  magnitudes 
show  semidiurnal  tidal  signals  and,  when  the  seagrass  bed  is  not  present, 
they  decrease  shoreward  due  to  increased  bottom  friction.  Current  magni¬ 
tudes  at  the  onshore  position  during  floods  are  slightly  smaller  than  ebbs 
because  flooding  tides  are  against  wave-induced  alongshore  currents  (to¬ 
ward  positive  y  direction).  However,  when  a  200-m-wide  bed  is  added 
(Eigure  24a),  current  magnitude  inside  the  bed  is  reduced  and  becomes 
smaller  than  either  offshore  or  onshore. 

In  the  time  series  of  five  other  variables,  tidal  signals  are  weak,  especially 
for  shallower  locations,  indicating  that  the  sediment  dynamics  in  the  sys¬ 
tem  are  dominated  by  waves.  Most  importantly,  averaged  suspended 
sediment  concentration,  skin  friction  shear  stress,  erosion  rate,  and  depo¬ 
sition  rate  are  lower  and  there  is  more  bottom  sediment  at  both  the  sea¬ 
grass  bed  and  onshore  positions  when  the  seagrass  bed  is  present.  This  re¬ 
sult  confirms  the  anticipation  of  lower  suspended  sediment  concentration 
from  Eigure  23  and  suggests  that  seagrass  beds  can  protect  bottom  sedi¬ 
ments  from  being  eroded  not  only  inside  the  bed  itself  but  also  the  area 
behind  it. 
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hours  hours 

Figure  24.  Time  series  of  six  variabies  for  the  configuration  of  200  m-wide  bed  (a)  and  no- 
seagrass  (b),  as  described  in  Scenario  4.  The  variabies  are  bottom  sediments  (kg/m^), 
suspended  sediment  concentration  (kg/m^),  skin  friction  shear  stress  (Pa),  magnitude  of 
current  (m/s),  erosion  rate  (kg/m^s),  and  deposition  rate  (kg/m^s)  from  the  top  to  bottom 
panei.  The  time  series  are  coiiected  aiongthe  central  transect,  as  shown  in  soiid  shore- 
normai  iine  in  Figure  24  (a).  The  soiid  iine  here  represents  average  of  each  variabie  from  the 
offshore  boundary  to  the  offshore  edge  of  the  bed  (450  m  from  offshore  boundary);  the 
dashed  iine  is  the  average  within  the  bed  (from  450  to  650  m);  the  dotted  iine  is  the  average 

over  the  rest  of  the  domain. 


Figure  25  shows  a  snapshot  of  distributions  of  suspended  (lower  panels) 
and  bottom  sediments  (upper  panels)  with  and  without  the  seagrass  bed  at 
maximum  flood.  For  the  no-seagrass  case  (Figure  25b),  suspended  sedi¬ 
ment  concentration  increases  shoreward  with  little  alongshore  variation, 
causing  bottom  sediments  to  decrease.  This  pattern  again  indicates  the 
dominance  of  wave-induced  erosion.  Adding  a  200-m-wide  seagrass  bed 
induces  both  alongshore  and  cross-shore  variations  of  suspended  and  bot¬ 
tom  sediment  distributions,  as  can  be  seen  in  Figure  25a.  Due  to  higher 
drag  of  the  bed,  tidal  currents  are  forced  to  flow  around  it,  resulting  in  a 
bulge  of  suspended  sediments  at  the  upstream  offshore  corner  of  the  bed. 
A  similar  pattern  is  observed  at  the  downstream  offshore  corner  when 
tides  change  direction.  In  general,  suspended  sediment  concentration 
within  the  bed  is  lower  than  that  either  onshore  or  offshore,  but  advection 
of  suspended  sediments  by  tidal  currents  can  locally  increase  the  concen¬ 
tration  within  the  bed.  As  for  bottom  sediments,  local  scouring  is  evident 
at  the  corners  of  the  bed  on  the  nearshore  side.  The  scouring  could  be  due 
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to  enhanced  tidal  flow  speed  between  the  shoreline  and  the  bed.  However, 
fleld  observations  could  not  be  found  to  support  such  scouring  in  the  lit¬ 
erature  and  suspect  that  this  effect  may  be  exaggerated  by  the  wall  bound¬ 
ary  condition  in  the  model.  Besides,  there  are  generally  more  bottom  sedi¬ 
ments  within  the  bed  than  on  either  side,  mostly  near  the  upper  and  lower 
edges.  The  sediment  trapping  is  due  to  import  of  higher  suspended  sedi¬ 
ment  concentration  by  tidal  currents  from  outside,  deposition  of  these 
sediments,  and  lower  wave-induced  erosion  rates  inside. 

Animating  the  model  results  confirms  that  sediment  trapping  appears  to 
occur  at  the  upstream  edge  on  each  half  tidal  cycle,  when  tidal  currents  are 
advecting  higher  suspended  sediment  concentrations  from  outside  into  the 
seagrass  bed. 

Implications  and  utility  of  modeling  studies 

Wave  attenuation  by  seagrass  beds 

Several  general  statements  can  be  made  from  the  results  of  the  model  sce¬ 
narios  with  a  flat  bottom  (Figures  20-22).  First,  larger  coverage  of  sea- 
grasses  in  the  direction  of  wave  propagation  results  in  higher  wave  at¬ 
tenuation  (indicated  by  percent  of  energy  flux  reduction)  and  thus  less 
energy  reaching  the  shoreline.  Second,  the  total  force  acting  on  the  bottom 
(indicated  by  percent  of  bottom  stress  reduction)  over  the  whole  domain 
decreases  as  the  seagrass  bed  is  moved  offshore.  Third,  relative  wave  at¬ 
tenuation  and  bottom  stress  reduction  are  greater  for  larger  incoming 
wave  heights. 

These  statements  are  generally  valid  as  long  as  there  is  significant  interac¬ 
tion  between  wave  orbital  velocities  and  the  seagrass  canopy.  This  qualifi¬ 
cation  may  be  interpreted  as  a  generalization  of  suggestions  by  Ward  et  al. 
(1984)  and  Koch  (2001).  They  pointed  out  that  wave  attenuation  should  be 
higher  when  seagrass  occupies  a  large  portion  of  the  water  column.  In¬ 
deed,  for  our  flat  bottom  cases,  seagrass  shoots  occupied  the  entire  water 
column  (June  case),  and  the  decay  of  orbital  velocity  with  depth  was  negli¬ 
gible  (at  i-m  depth,  4  sec  waves  are  close  to  shallow- water  waves).  In  con¬ 
trast,  in  our  sloping  bottom  cases  i-m  seagrass  shoots  only  occupied  part 
of  the  water  column  in  the  deepest  region  (2.5-m  depth),  and  orbital  veloc¬ 
ity  decayed  at  least  25  percent.  This  is  why  the  differences  in  wave  height 
between  no-seagrass  and  full-coverage  cases  (Figure  8a)  are  relatively 
small  in  the  deepest  region  but  increase  in  shallower  regions.  Although  a 
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Figure  25.  Snapshot  (top  view  of  modei  domain)  of  distribution  of  bottom  sediments  (upper 
panei)  and  suspended  sediment  concentration  (bottom  panei)  for  the  configuration  of 
200-m-wide  bed  (a)  and  no-seagrass  (b).  The  M2  tide  is  forced  in  the  aiongshore  direction, 
whiie  0.1  m,  4  sec  waves  propagate  from  offshore  boundary  with  10  deg  incident  angie 
(counterciockwise  from  shore-normai  direction).  The  current  direction  and  magnitude  are 
indicated  by  vectors,  and  bottom  sediment  and  suspended  sediment  concentration  are 
shown  in  contour.  The  soiid  iine  in  (a)  represents  the  centrai  transect  aiong  which  model 
outputs  in  Figure  7  and  8  are  coiiected.  Itshouid  be  noted  that  the  aiongshore  domain 
(5,400  m)  is  set  smaiier  than  the  domain  of  SHORECiRC  /  REFDIF,  and  a  iooping 
boundary  condition  is  appiied. 
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flume  study  (Fonseca  and  Cahalan  1992)  and  field  observations  (Koch 
1996)  indirectly  support  this  hypothesis,  wave  attenuation  has  also  been 
observed  when  seagrass  only  occupies  a  small  portion  of  the  water  column 
(Granata  et  al.  2001;  depth  of  15  m).  Systematic  observations  on  the  effects 
of  seagrass  bed  geometry  on  waves  with  different  wave  heights  and  periods 
are  needed  to  verify  model  predictions  and  to  better  understand  these 
processes. 

Wave  attenuation  by  seagrass  may  have  implications  for  shoreline  protec¬ 
tion  (van  Katwijk  2000;  Koch  et  al.  2006).  Using  observed  seagrass  pa¬ 
rameters  in  June,  our  model  results  showed  significant  reductions  in  wave 
energy  flux  at  the  shoreline  for  both  flat  and  sloping  bottom  cases  (Figure 
20,  21,  and  23).  However,  seagrasses  vary  seasonally  in  temperate  envi¬ 
ronments,  whereas  shoreline  erosion  is  usually  associated  with  large  wave 
events  that  occur  episodically  (Wilcock  et  al.  1998)  over  annual  or  decadal 
time  scales  (Kamphuis  1987).  Hence,  the  timing  between  wave  events  and 
seagrass  growth  influences  the  potential  for  seagrass  beds  to  protect 
shorelines.  Without  knowing  more  about  this  timing,  it  is  difficult  to 
evaluate  the  net  influence  of  seagrass  on  shoreline  protection  based  on  the 
results  presented  here.  Other  factors  such  as  the  spectral  or  directional 
distributions  of  wave  energy  may  also  need  to  be  considered  in  order  to 
better  address  this  question.  REF/DIF  is  capable  of  modeling  spectral 
wave  forcing  as  well  as  multiple  wave  incident  angles  (Kirby  and  Tuba 
Ozkan  1994),  closer  to  what  is  observed  in  the  field.  Those  issues  will  be 
addressed  in  the  future. 

Sediment  dynamics 

Model  results  have  two  main  implications  for  sediment  dynamics.  First, 
sedimentary  processes  are  altered  within  the  seagrass  bed  and  probably 
behind  it.  Results  from  Scenario  3  (Figure  22)  show  that,  in  wave- 
dominated  environments,  the  total  force  acting  on  bottom  sediments  de¬ 
creases  as  the  seagrass  bed  is  moved  offshore  due  to  increases  in  the  af¬ 
fected  area  behind  the  bed.  This  suggests  that  seagrass  beds  may  stabilize 
bottom  sediments  in  the  zone  between  the  bed  and  shoreline,  which  is 
consistent  with  Hine  et  al.  (i987)’s  observation  that  disappearance  of  a 
seagrass  community  allowed  rapid  onshore  and  longshore  sand  transport 
in  the  nearshore  zone.  Similarly,  comparison  between  cases  with  and 
without  a  200-m-wide  bed  (Figure  24)  shows  lower  skin  friction  shear 
stresses,  lower  erosion  rates,  and  higher  levels  of  bottom  sediments  at  lo¬ 
cations  within  and  behind  the  bed.  Within  the  bed,  the  result  is  consistent 
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with  Lopez  and  Garcia  (i998)’s  findings.  Applying  a  numerical  model  of 
open  channel  flow  with  simulated  vegetation,  they  showed  reduced  shear 
stress  and  consequently  lower  suspended  sediment  transport  (partly  due 
to  lower  suspended  sediment  concentration)  in  the  vegetated  area.  Re¬ 
duced  erosion  rate  as  well  as  bottom  sediment  retention  are  also  sup¬ 
ported  by  field  observations  (Gacia  and  Duarte  2001).  Gacia  and  Duarte 
found  that  the  presence  of  Posidonia  oceanica  enhances  sediment  stability 
by  preventing  resuspension.  Behind  the  bed,  however,  there  is  no  quanti¬ 
tative  evidence  to  support  the  model-predicted  reduction  in  skin  friction 
shear  stress  and  erosion  and  sediment  retention.  Further  studies  are  re¬ 
quired  for  verification. 

Second,  sediment  trapping  in  the  seagrass  bed  requires  horizontal  trans¬ 
port  to  bring  suspended  sediment  from  outside  into  the  bed.  The  concept 
of  a  seagrass  bed  as  a  depositional  environment  has  been  suggested  by 
several  authors  (e.g.,  Grady  1981;  Ward  et  al.  1984;  Almasi  et  al.  1987), 
and  the  proposed  mechanism  for  this  accumulation  may  be  summarized  as 
reduced  shear  stress  due  to  loss  of  momentum  in  a  seagrass  bed  leading  to 
reduction  in  resuspension  and  thus  increased  sediment  accumulation 
(Koch  et  al.  2006).  The  connection  between  lower  momentum  and  re¬ 
duced  resuspension  (lower  erosion  rate)  is  supported  by  our  results.  How¬ 
ever,  the  results  suggest  that  linkages  from  reduced  resuspension  to  in¬ 
creased  accumulation  are  not  trivial  and  may  not  occur  everywhere  within 
the  bed.  Sediment  accumulation  occurs  when  the  suspended  sediment 
concentration  is  high  enough  that  the  deposition  rate  exceeds  the  erosion 
rate.  Sediment  accumulation  at  the  upper  and  lower  edges  of  the  bed  in 
Figure  25a  illustrates  this  point.  At  these  two  edges,  accumulation  occurs 
when  high  concentrations  of  suspended  sediment  from  outside  are  trans¬ 
ported  into  the  bed,  where  reduced  shear  stresses  allow  them  to  deposit. 
The  amount  of  accumulated  sediment  then  gradually  decreases  with  dis¬ 
tance  into  the  bed  until  the  sediment  source  from  outside  is  used  up.  Be¬ 
yond  this  location,  sediments  that  were  originally  there  remain,  but  there 
is  no  new  accumulation. 

In  short,  the  model  results  clearly  demonstrate  that  sediment  accumula¬ 
tion  requires  both  sediment  sources  (outside  the  bed)  and  a  transport 
mechanism  (tidal  currents),  both  of  which  may  vary  spatially  within  the 
bed.  The  reduction  in  suspended  sediment  transport  capacity  (concentra¬ 
tion  multiplied  by  streamwise  velocity)  in  a  vegetated  area  reported  by 
Lopez  and  Garcia  (1998)  indirectly  supports  accumulation  at  the  bed  edge. 
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Nevertheless,  direct  observations  of  spatial  patterns  of  accumulation 
within  seagrass  beds  are  few,  and  most  of  them  focus  on  sediment  grain 
size  distributions  (e.g.,  Scoffin  1970;  Wanless  1981;  Granata  et  al.  2001). 
More  observations  that  can  resolve  spatial  patterns  of  erosion/deposition 
are  needed  to  move  our  understanding  forward.  It  should  also  be  noted 
that  the  spatial  distribution  of  bottom  sediment  presented  here  may  not 
match  field  observations  exactly  because  these  model  scenarios  do  not  ac¬ 
count  for  limited  supplies  of  surficial  sediments  or  mixed  sediment  grain 
sizes. 

Utility  of  modeling  studies 

Teeter  et  al.  (2001)  summarized  the  complexities  of  modeling  interactions 
between  waves,  flow,  and  sediment  transport  in  seagrass  beds  quite 
clearly,  and  pointed  out  some  of  the  limitations  in  our  current  capabilities. 
The  model  presented  here  underscores  the  importance  of  adequately  rep¬ 
resenting  these  complex  interactions.  For  example,  it  is  quite  clear  from 
the  model  runs  that  wave  forcing  alone  will  not  result  in  enhanced  deposi¬ 
tion  within  a  seagrass  bed.  Wave  forcing  can  create  an  erosional  environ¬ 
ment  outside  and  a  depositional  environment  inside,  but  without  currents 
forced  by  winds  or  tides,  there  will  be  little  transport  of  suspended  sedi¬ 
ment  from  outside  to  inside.  It  is  important  to  point  out  that  our  model  in 
fact  does  not  represent  all  of  the  potential  interactions  identified  by  Teeter 
et  al.  (2001),  ignoring  direct  wind  forcing,  multiple  sediment  grain  sizes, 
and  large  tidal  height  variability;  and  neither  the  framework  presented  by 
Teeter  et  al.  nor  the  present  model  consider  ecosystem  effects  (such  as 
feedbacks  to  seagrass  growth). 

Given  the  complexities  of  these  interactions,  the  uncertainties  associated 
with  attempting  to  model  some  of  them,  and  the  daunting  task  of  develop¬ 
ing  and  successfully  implementing  the  equally  complex  computer  codes 
required,  it  might  be  asked  if  the  effort  is  justified.  The  answer  from  our 
point  of  view  is  a  resounding  “yes,”  but  not  in  isolation.  We  hope  that  some 
of  the  example  scenarios  presented  here  have  illustrated  the  utility  of  the 
modeling  approach,  especially  in  combination  with  sparse  data.  The  pri¬ 
mary  advantages  of  modeling  are: 

a.  To  generate  testable  hypotheses  or  to  help  explain  observations.  For 
example,  the  model  that  was  presented  predicts  that  sediment 
deposition  in  seagrass  beds  should  be  most  pronounced  just  inside 
the  bed  on  the  upstream  side,  and  that  internal  gradients  in  deposi- 
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tion  will  be  quite  large.  Is  this  observable?  The  model  also  predicts 
that  wave  attenuation  expressed  as  a  percentage  of  the  incoming 
wave  energy  will  be  greater  for  higher  waves,  which  is  apparent 
qualitatively  in  data  presented  elsewhere  in  this  report. 

b.  To  aid  in  the  design  of  observational  programs  or  restoration  ef¬ 
forts.  From  an  experimental  design  perspective,  given  reasonable 
estimates  of  tidal  currents,  waves,  seagrass  densities,  and  bed  sizes, 
a  model  might  be  used  to  anticipate  where  sparse  instrumentation 
could  be  put  to  best  use,  and  what  kind  of  measurements  might  be 
expected  to  yield  the  most  useful  observations.  From  a  restoration 
point  of  view,  if  the  goal  is  to  use  seagrass  beds  to  reduce  nearshore 
turbidity  due  to  resuspension.  Figure  22  indicates  that  a  bed  placed 
further  offshore  will  have  a  greater  effect  on  total  bottom  stress  re¬ 
duction  (hence,  total  resuspension)  -  that  is,  unless  the  bed  is  also 
so  deep  that  it  falls  below  wave  base  and  no  longer  affects  the 
waves.  However,  the  bed  distance  offshore  has  relatively  little  pre¬ 
dicted  impact  on  reduction  in  wave  power  striking  the  shoreline 
(not  shown). 

c.  To  help  interpolate  sparse  data  in  space  and  time.  We  can  never 
measure  enough  to  cover  all  scales  of  interest  and  all  times,  but  if 
we  can  measure  enough  to  build  and  parameterize  a  reliable  model, 
we  can  fill  in  the  gaps. 

For  these  reasons,  we  advocate  the  development  and  use  of  numerical 
models  in  conjunction  with  careful  field  and  laboratory  measurements  as 
the  best  possible  combination  of  tools  to  address  both  research  and  resto¬ 
ration  goals. 
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7  Summary  and  Conclusions 

While  extensive  work  has  been  done  on  how  currents  move  through  sea- 
grass  beds,  how  currents  are  affected  by  seagrasses  and  how  seagrasses  are 
affected  by  currents,  our  understanding  about  how  waves  propagate 
through  seagrass  beds,  are  affected  by  seagrasses  and  affect  seagrass  biol¬ 
ogy  is  still  in  its  infancy.  It  has  been  just  recently  that  more  attention  has 
been  paid  to  this  field  of  study.  By  using  experimental  designs  and  instru¬ 
mentation  that  allows  for  easy  comparison  between  studies,  it  will  be  pos¬ 
sible  to  gather  the  information  so  much  needed  to  understand  waves  in 
vegetated  systems.  The  continuous  development  of  smaller  electronics 
now  also  allows  for  the  deployment  of  wave  sensors  in  relatively  shallow 
waters.  Modeling  of  waves  in  seagrass  colonized  areas  has  also  seen  a  tre¬ 
mendous  leap  in  recent  years  but  many  questions  still  remain  to  be  ex¬ 
plored.  The  continuous  and  careful  collection  of  field  and  lab  data  on 
waves  in  seagrass  beds  of  different  species,  at  different  depths  and  at  dif¬ 
ferent  times  of  the  year  will  allow  existing  models  to  be  refined  and  new 
models  to  be  developed.  In  summary,  wave-related  processes  in  seagrass 
systems  is  a  fascinating  field  in  which  extensive  contributions  can  still  be 
made. 
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Appendix;  Glossary  of  Terms 


Acclimation 

Process  by  which  plants  and  animals  physiologically  adapt  to  envi¬ 
ronmental  conditions.  It  differs  from  ‘adaptation’  in  that  it  does  not 
involve  genetics. 

Acoustic  Doppler 
Current  Profiler 

An  instrument  that  uses  the  Doppler  shift  of  reflected  sound  from 
projected  acoustic  beams  to  measure  profiles  of  velocity,  usually 
from  Just  above  the  bottom  to  Just  below  the  surface. 

Advection 

The  transport  of  a  dissolved  or  suspended  quantity  by  flowing  water 
(here  it  refers  to  horizontal  transport  of  sediments  by  currents). 

Algae 

Aquatic  plants  (freshwater  to  marine)  that  reproduce  via  spores.  Al¬ 
gae  can  be  composed  of  a  single  cell  or  reach  25  m  in  length  (kelp). 
Also  see  ‘Macroalgae,’  Phytoplankton’  and  ‘seaweed.’ 

Bathymetry 

The  spatial  distribution  of  bottom  depth  in  bodies  of  water. 

Biomass 

Weight  of  biological  material  (plants  or  animals)  expressed  as  a 
function  of  area  or  volume. 

Bottom  roughness 
coefficient 

Within  the  bottom  boundary  layer,  the  value  of  the  0  flow  speed  in¬ 
tercept  height  on  a  plot  of  flow  speed  against  the  logarithm  of  height 
above  bottom.  In  rough  turbulent  flow,  the  bottom  roughness  coeffi¬ 
cient  is  equal  to  the  physical  bottom  roughness  divided  by  30.  A 
larger  roughness  coefficient  means  higher  drag  on  flow,  given  the 
same  flow  speed. 

Bottom  shear  stress 

The  force  per  unit  area  acting  on  the  bottom  due  to  strong  vertical 
gradients  (shear)  in  the  fluid  immediately  above  the  bottom.  Bottom 
shear  stress  is  often  computed  as  a  drag  coefficient  multiplied  by 
velocity  (at  a  defined  height)  squared. 

Bottom  sheltering 

Protection  of  the  bottom  from  fluid  forces,  caused  by  the  form  drag 
of  obstacles  such  as  seagrasses  and/or  the  deflection  of  flow  away 
from  the  bottom. 

Boundary  layer 

The  layer  of  reduced  velocity  that  is  immediately  adjacent  to  the  sur¬ 
face  of  a  solid  past  which  the  water  is  flowing. 

Breaking  zone 

The  nearshore  zone  in  which  wave  heights  decrease  rapidly  and  tur¬ 
bulence  levels  increase  markedly  due  to  wave  breaking. 

Breakwater 

A  barrier  that  protects  the  shoreline  from  the  full  impact  of  waves 

Canopy  height 

Height  above  the  sediment  to  which  most  (2/3)  of  the  seagrass 
leaves  grow  (see  Figure  2). 

Chop  or  wind-sea 

Short,  steep  water  waves  of  higher  frequency  generated  locally  by 
the  wind,  with  a  broad  distribution  of  wavelengths/periods. 

Cohesive  sediment 

Sediments  composed  of  clay,  silt,  and/or  organic  matter  for  which 
the  dominant  stabilizing  force  is  the  tendency  of  particles  to  stick 
together,  rather  than  the  weight  of  the  particles. 
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Consolidation 

The  process  by  which  a  sediment  layer  decreases  in  height  due  to 
the  expulsion  of  water.  Consolidation  often  occurs  when  the  weight 
of  a  sediment  layer  squeezes  the  water  out  of  the  space  between 
sediment  particles. 

Crest 

The  highest  point  of  a  wave  cycle. 

Curvilinear 

Describes  a  coordinate  system  that  bends  to  follow  the  shape  of  the 
bottom  or  coastline,  as  opposed  to  a  Cartesian  coordinate  system 
with  only  rectulangular  cells. 

Diatom 

Single  celled  phytoplankton  with  an  external  skeleton  of  silica. 

Diffusive  boundary 
layer 

Sublayer  of  the  boundary  layer  where  the  flux  of  molecules  occurs 
via  diffusion. 

Directional  wave 
spectrum 

The  distribution  of  wave  energy  as  a  function  of  both  wave  frequency 
(period'^)  and  wave  direction  (usually  the  direction  from  which  the 
waves  are  coming). 

DIWASP 

A  Directional  WAve  SPectra  toolbox  for  MATLAB. 

Drag  coefficient 

The  ratio  between  fluid  drag  and  the  product  of  fluid  density  and 
velocity  squared,  for  flow  over  a  boundary  or  around  an  obstacle;  a 
higher  drag  coefficient  means  greater  drag  for  the  same  velocity. 

Duration 

The  amount  of  time  for  which  a  given  wind  has  been  blowing. 

Ecosystem  engineer 

An  organism  that  has  the  capacity  to  alter  the  environment  it  inhab¬ 
its. 

Ecosystem  service 

Goods  provided  by  natural  communities  (e.g.  the  capacity  of  sea- 
grasses  to  remove  nutrients  from  the  water  column). 

Ebb 

Outflow  of  the  tide. 

Electromagnetic 
velocity  sensor 

An  instrument  that  measures  components  of  water  velocity  at  a  fixed 
point  in  space  by  sensing  flow-induced  changes  in  an  emitted  elec¬ 
tromagnetic  field. 

Energy  dissipation 

Loss  of  energy  over  time  due  to  breaking  or  friction. 

Epiphyte 

Plants  that  grow  on  the  surface  of  other  plants.  In  the  case  of  sea- 
grasses,  algae  that  grow  on  seagrass  leaves  and  host  a  community 
of  associated  bacteria  and  grazers. 

Erosion  rate 

The  rate  at  which  sediment  is  eroded  (or  resuspended)  from  the  bot¬ 
tom. 

Euphotic  zone 

Surface  waters  that  receive  sufficient  light  to  support  photosynthesis 
(1  percent  of  surface  light). 

Fetch 

The  length  of  water  over  which  a  given  wind  has  blown. 

Fine  particles 

Particles  with  sediment  grain  size  smaller  than  63  pm  (silt  and  clay). 

Flood 

Inflow  of  the  tide. 

Flow  speed 

The  rate  of  motion  of  a  fluid. 

Flux  boundary 
condition 

A  type  of  model  boundary  condition  in  which  the  time-dependent 
fluid  volume  flux  (flow  speed  times  cross-sectional  area)  across  the 
boundary  is  specified. 
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Form  drag 

Frictional  force  that  arises  from  the  form  of  an  object,  in  addition  to 
the  drag  of  its  surface.  Form  drag  is  primarily  due  to  flow  separation 
and  turbulence  behind  the  object. 

Free-stream  velocity 

A  reference  velocity  away  from  the  region  of  boundary  or  object  in¬ 
fluence.  Flere  it  refers  to  the  axial  velocity  in  a  flume  measured  up¬ 
stream  of  the  leading  edge  of  a  seagrass  bed  (so  not  influenced  by 
the  seagrass). 

Frictional  damping 

Momentum  or  energy  loss  due  to  friction. 

Frictional  drag 

Force  that  opposes  fluid  motion  due  to  flow  over  a  boundary  or 
around  an  obstacle. 

Geomorphology 

Description  of  the  nature  and  history  of  the  landforms  and  the  proc¬ 
esses  that  create  them. 

Group  velocity 

The  rate  at  which  wave  energy  propagates  through  space. 

Flabitat 

The  place  inhabited  by  a  plant  or  animal. 

Fleaviside  step 
function 

In  the  present  context,  a  function  defined  to  be  0  for  all  arguments 
less  than  0,  and  1  for  all  arguments  greater  than  or  equal  to  0. 

Infragravity  wave 

Long  waves  with  periods  of  30  sec  to  several  minutes,  usually  gen¬ 
erated  by  interactions  between  wave  refraction/breaking  and  near¬ 
shore  bathymetry,  propagating  alongshore  rather  than  across-shore. 

Kelp 

A  brown  macroalga  found  in  relatively  cold  waters  that  has  a  distinct 
morphology  with  leaves  (fronds)  that  can  reach  75  ft  in  length. 

Kelvin-  Flelmholtz 
instabilities 

Instability  that  occurs  when  shear  is  present  within  a  continuous 
fluid  or  when  there  is  sufficient  velocity  difference  across  the  inter¬ 
face  between  two  fluids  (e.g.,  between  the  seagrass  canopy  and  the 
water  above  it). 

r 

1 

Source:  httD://www.answers.com/toDic/wavecloudsduval-iDf‘. 

Kinematic  viscosity 

Measure  of  the  resistance  of  a  fluid  to  shear  deformation  under  an 
applied  force,  normalized  by  fluid  density. 

Light-fleck 

Pulse  or  speck  of  light.  In  seagrass  habitats,  light-flecks  are  the  re¬ 
sult  of  light  focusing  by  the  crest  of  waves  in  shallow  wate.r 

Linear  or  airy  wave 
approximation 

An  approximation  assuming  that  wave  amplitude  is  small,  which 
greatly  simplifies  the  equations  of  wave  motion. 
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Macroalga 

Large  (visible  to  the  naked  eye)  aquatic  plant  (freshwater  to  marine) 
that  reproduces  via  spores  and  has  a  relatively  simple  morphology: 
no  roots,  rhizomes  or  flowers  (also  referred  to  as  ‘seaweed’).  Found 
on  rocky  substrates. 

MATLAB 

A  high-level  numerical  computing  environment  and  programming 
language  created  by  The  MathWorks,  Inc.,  Natick,  Massachusetts. 

Microalga 

Small  (unicellular  or  clusters  of  cells)  aquatic  plant  (freshwater  to 
marine). 

Model  domain 

The  equivalent  physical  space  of  a  model  simulation.  Here,  the 
model  domain  is  a  rectangular  area  720m  x  5400  m. 

Monochromatic 

Having  a  single  wavelength  or  wave  period. 

Nearshore 

The  zone  that  extends  from  the  shoreline  to  a  position  marking  the 
start  of  the  offshore  zone,  typically  at  water  depths  of  approximately 
3-5  m. 

Organic  content 

Amount  of  organic  matter  found  in  sediments.  Often  determined  as 
the  fraction  of  the  sediment  that  is  combustible  at  450  °C. 

Oscillatory  flow 

The  back-and-forth  motion  of  a  fluid  under  waves. 

Period 

The  time  for  one  complete  cycle  of  a  wave  to  pass  by  a  fixed  point  in 
space. 

Phase  speed 

The  rate  at  which  a  wave  pattern  propagates  through  space. 

Photosynthesis 

Process  by  which  plants  transform  carbon  dioxide,  solar  energy  and 
water  into  chemical  energy  needed  for  growth. 

Phytoplankton 

Microscopic  algae  found  in  the  euphotic  zone. 

Pressure 

The  force  per  unit  area  applied  on  a  surface  in  a  direction  perpen¬ 
dicular  to  that  surface. 

Productivity 

The  rate  at  which  a  given  quantity  of  organic  material  is  produced  by 
organisms  such  as  seagrasses. 

Recruitment 

The  addition  of  new  individuals  to  a  population. 

REF/DIF 

A  model  for  ocean  surface  wave  propagation  that  accounts  for  wave 
bending  when  waves  travel  through  an  area  with  changing  depths 
(refraction)  and  wave  scattering  behind  obstacles  (diffraction).  It  was 
developed  at  the  University  of  Delaware. 

Resuspension 

The  process  whereby  previously  deposited  sediment  particles  are 
lifted  back  into  the  water  column  by  fluid  (or  biogenic)  forces. 

Reynolds  number 

The  ratio  of  inertial  to  viscous  force  in  a  fluid  flow;  low  Reynolds 
numbers  correspond  to  laminar  flow  (smooth,  steady  streamlines) 
while  high  Reynolds  numbers  correspond  to  turbulent  flow  (smooth 
average  streamlines,  but  chaotic  instantaneous  flow). 

Rhizome 

Horizontal  structure  of  seagrasses  that  connects  shoots  and  roots. 

Rough  turbulent  flow 

A  fully  turbulent  flow,  such  that  increasing  the  Reynolds  number 
does  not  change  the  qualitative  behavior  of  the  turbulence.  This 
means  that  molecular  viscosity  no  longer  plays  a  role;  for  example,  in 
rough  turbulent  boundary  layer  flow  there  is  no  viscous  sublayer  just 
above  the  boundary. 
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Roughness 

Measurement  of  the  small-scale  variations  in  the  height  of  a  surface. 
Higher  bottom  roughness  often  means  larger  fluid  drag  and  total 
bottom  stress. 

Sandbar 

Submerged  (or  partly  submerged)  ridge  of  sand  along  a  shore. 

Sand  wave 

Large  regular  alterations  of  sand  ridges  (crests  and  troughs)  which 
move  when  current  velocities  exceed  47  to  60  cm/sec  causing  a 
succession  of  depositional  and  erosional  events  at  a  fixed  point. 

SAV 

See  ‘submersed  aquatic  vegetation.’ 

Scour 

Preferential  local  erosion  of  bottom  sediment  due  to  focused  cur¬ 
rents,  waves,  or  turbulence. 

Sea-surface  elevation 

The  height  of  the  sea-surface  with  respect  to  a  fixed  reference  height 
such  as  mean  sea-level. 

Sea-surface  slope 

The  inclination  of  the  sea-surface  from  a  perfectly  horizontal  line. 

Sea-surface  vertical 
acceleration 

Rate  of  change  of  the  sea-surface  vertical  velocity. 

Sea-surface  vertical 
velocity 

Velocity  component  of  the  sea-surface  in  the  vertical  direction. 

Seagrass 

Flowering  marine  plant. 

Seaweed 

See  ‘macroalgae.’ 

Sediment  deposition 

The  process  whereby  sediments  in  the  water  column  sink  to  the  bot¬ 
tom  and  attach  to  it. 

Sediment  dynamics 

The  study  of  the  effects  of  physical  (and  sometimes  biological)  forces 
on  the  motion  of  sediments. 

Sediment  grain  size 

Characteristic  physical  dimension  of  sediment  particles.  For  exam¬ 
ple,  a  grain  size  for  fine  sand  might  be  0.2  mm. 

Sediment  transport 

The  motion  of  sediment  in  a  flowing  fluid,  often  separated  into  sus¬ 
pended  transport  and  bedload  transport.  Suspended  transport  de¬ 
scribes  the  3-D  movement  of  sediments  that  stay  in  the  water  col¬ 
umn;  Bed-load  transport  describes  the  rolling  and  hopping  of 
sediment  particles  across  the  bottom. 

Self-shading 

Process  of  one  part  of  the  plant  shading  the  same  plant  (e.g.  when  a 
seagrass  leaf  bends  over  it  shades  the  leaf  underneath). 

Semidiurnal 

Occurring  approximately  once  every  12  hr. 

Separation 

When  a  flow  streamline  detaches  from  the  surface  of  an  object, 
such  that  the  mean  flow  beyond  that  point  is  zero  or  recirculating. 

Flow  beyond  the  point  of  separation  is  dominated  by  turbulence,  so 
early  flow  separation  usually  means  higher  drag. 

Shallow  water  wave 

Waves  for  which  the  local  water  depth  is  less  than  5%  of  the  wave¬ 
length. 

Shear  velocity 

The  main  scale  velocity  for  boundary  layer  turbulence,  defined  as  the 
square  root  of  bottom  shear  stress  divided  by  fluid  density. 

Shoaling 

Changes  in  wave  height  and  wavelength  as  waves  propagate  into 
shallower  water. 
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Shoot 

Cluster  of  leaves  (see  Figure  1). 

SHORECIRC 

A  numerical  nearshore  circulation  model  that  accounts  for  currents 
induced  by  wave  breaking,  as  well  as  other  more  conventional 
forces.  SHORECIRC  was  developed  at  the  University  of  Delaware. 

Sill 

Submersed  breakwater. 

Skin  friction 

Frictional  force  that  arises  from  shear  stresses  generated  at  the  sur¬ 
face  of  an  object  or  at  the  bottom  of  a  water  column.  Skin  friction  is 
the  force  that  most  directly  dislodges  sediment  particles  or  controls 
interfacial  fluxes. 

Steady  current 

Flow  speed  that  doesn’t  change  with  time. 

Submersed  aquatic 
vegetation  (SAV) 

Flowering  plants  growing  submersed  in  habitats  covering  a  wide 
range  of  salinities,  from  freshwater  to  marine.  Note  that  seagrasses 
can  also  be  referred  to  as  SAV. 

Surface  gravity  wave 

Waves  traveling  on  the  free  surface  of  a  fluid  medium  and  having  a 
restoring  force  of  gravity  or  buoyancy. 

Suspended  sediment 
concentration 

The  mass  of  suspended  sediment  in  a  unit  volume  of  water. 

Swells 

Long  period,  long  crested  surface  waves  generated  far  away  and 
propagating  into  a  region  of  interest. 

Tidal  current 

Oscillatory  velocity  induced  by  the  surface  tide,  with  flood  corre¬ 
sponding  to  the  incoming  tide  and  ebb  corresponding  to  the  outgo¬ 
ing  tide. 

Transect 

A  reference  line  drawn  through  a  domain  of  interest.  Here  it  refers  to 
a  slice  that  is  perpendicular  to  the  shoreline. 

Trough 

The  lowest  point  of  a  wave  cycle. 

Turbulent  wake 
interference 

Occurs  when  the  turbulent  wake  of  an  upstream  object  encounters  a 
downstream  object. 

Vector  addition 

Summation  of  two  vectors  by  the  addition  of  their  directional  com¬ 
ponents;  for  example,  adding  two  east  components  of  velocity  and 
two  north  components  of  velocity  to  calculate  the  total  velocity  vec¬ 
tor. 

Von  Karman  constant 

The  constant  of  proportionality  between  shear  velocity  and  the  prod¬ 
uct  of  height  and  velocity  shear  at  any  height  z  in  the  constant  stress 
region  of  a  turbulent  boundary  layer.  It  was  determined  experimen¬ 
tally  to  have  a  value  of  0.4. 

Wave 

A  disturbance  that  propagates  through  space,  transmitting  energy 
but  relatively  little  mass. 

Wave  energy 
attenuation 

Decrease  in  wave  energy,  as  the  wave  propagates  due  to  frictional 
forces  or  spreading  of  wave  energy,  or  with  depth  below  the  surface 
for  deep  water  or  intermediate  waves. 

Wave  energy  flux  or 
wave  power 

The  product  of  wave  energy  and  wave  group  velocity.  When  wave 
power  is  dissipated  by  bottom  friction  or  breaking,  the  wave  power 
lost  is  equal  to  the  work  done  by  the  waves  on  the  bottom  and/or 
shoreline. 
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